THE BIANCHI IX ATTRACTOR 



HANS RINGSTROM 



- z - Abstract. We consider the asymptotic behaviour of spatially homogeneous 

, spacetimes of Bianchi type IX close to the singularity (we also consider some of 

^—■^ i the other Bianchi types, e. g. Bianchi VIII in the stiff fluid case). The matter 

^ ' content is assumed to be an orthogonal perfect fluid with linear equation of 

(N ■ state and zero cosmological constant. In terms of the variables of Wainwright 

f*H , and Hsu, we have the following results. In the stiff fluid case, the solution 

■ converges to a point for all the Bianchi class A types. For the other matter 

models we consider, the Bianchi IX solutions generically converge to an at- 
tractor consisting of the closure of the vacuum type II orbits. Furthermore, 
we observe that for all the Bianchi class A spacetimes, except those of vacuum 
Taub type, a curvature invariant is unbounded in the incomplete directions of 
inextendible causal geodesies. 
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\^ . 1. Introduction 

o . 

, The last few decades, the Bianchi IX spacetimes have received considerable atten- 

' tion, see for instance and references therein. Agreement has been 

, reached, at least concerning some aspects of the asymptotic behaviour as one ap- 

O^' proaches a singularity, but the basis for the consensus has mainly consisted of nu- 

^ , merical studies and heuristic arguments. The objective of this article is to provide 

mathematical proofs for some aspects of the 'accepted' picture. The main result of 
this paper was for example conjectured in jisj p. 146-147, partly on the basis of a 
numerical analysis. 



X 



' Why Bianchi IX? One reason is the fact that this class contains the Taub-NUT 



spacetimes. These spacetimes are vacuum maximal globally hyperbolic spacetimes 
that are causally geodesically incomplete both to the future and to the past, see 
|6) and 10. However, as one approaches a singularity, in the sense of causal ge- 
odesic incompleteness, the curvature remains bounded. In fact, one can extend 
the spacetime beyond the singularities in inequivalent ways, see It is natural 
to conjecture that the behaviour exhibited by the Taub-NUT spacetimes is non- 
generic, and it is interesting to try to prove that the behaviour is non-generic in the 
Bianchi IX class. In fact we prove that all Bianchi IX initial data considered in this 
paper other than Taub-NUT yield inextendible globally hyperbolic developments 
such that the curvature becomes unbounded as one approaches a singularity. This 
result is in fact more of an observation, since the corresponding result is known in 



the vacuum case, see [16 , and curvature blow up is easy to prove in the non- vacuum 
cases we consider. 

Another reason for studying the Bianchi IX spacetimes is the BKL conjecture, 
see 1^. According to this conjecture, the 'local' approach to the singularity of 
a general solution should exhibit oscillatory behaviour. The prototypes for this 
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behaviour among the spatiaUy homogeneous spacetimes are the Bianchi VIII and 
IX classes. Furthermore the matter is conjectured to become unimportant as one 
approaches a singularity, with some exceptions, for example the stiff fluid case. 
We refer to [Q for arguments supporting the BKL conjecture and to for an 
overview of conjectures and results under symmetry assumptions of varying degree. 
In this paper we prove, under certain restrictions on the allowed matter models, 
that generic Bianchi IX solutions exhibit oscillatory behaviour and that the matter 
becomes unimportant as one approaches a singularity. What is meant by the latter 
statement will be made precise below. If the matter model is a stiff fluid the matter 
will be important, and in that case we prove that the behaviour is quiescent. This 
should be compared with |^ concerning the structure of singularities of analytic 
solutions to Einstein's equations coupled to a scalar field or stiff fluid. In that 
paper, Andersson and Rendall prove that given a certain kind of solution to the so 
called velocity dominated system, there is a unique solution of Einstein's equations 
coupled to a stiff fluid approaching the velocity dominated solution asymptotically. 
One can then ask the question whether it is natural to assume that a solution has 
the asymptotics they prescribe. In Section we show that all Bianchi VIII and 
IX stiff fluid solutions exhibit such asymptotic behaviour. 

The results presented in this paper can be divided into two parts. The first part 
consists of statements about developments of orthogonal perfect fluid data of class 
A. We clarify below what we mean by this. The results concern curvature blow up 
and inextendibility of developments. The second part consists of results expressed 
in terms of the variables of Wainwright and Hsu. These variables describe the 
spacetime close to the singularity, and we prove that Bianchi IX solutions generically 
converge to a set on which the flow of the equation coincides with the Kasner map. 

We consider spatially homogeneous Lorentz manifolds (M, g) with a perfect fluid 
source. The stress energy tensor is thus given by 

(1.1) Tab ^ fJ-UaUb + p{gab + UaUb), 

where m is a unit timelike vectorfield, the 4-velocity of the fluid. We assume that p 
and /i satisfy a linear equation of state 

(1.2) p={^~l)fi, 

where we in this paper restrict our attention to 2/3 < 7 < 2. We will also assume 
that u is perpendicular to the hypersurfaces of homogeneity. Einstein's equations 
can be written 

(1-3) Rab — -^Rfjab — Tab, 

where Rab and R are the Ricci and scalar curvature of (M, g). In order to formulate 
an initial value problem in this setting, consider a spacelike submanifold (M, g) of 
{M,g), orthogonal to u. Let Cq, a = 0, ..,3 be a local frame with cq = u and e^, 
i = 1,2,3 tangent to M and let fey be the second fundamental form of {M,g). 
Then g and k must satisfy the equations 

Rg - + (tr gkf = 2i?oo + R 

and 

Vjtrgfc - V^/cy = Roi, 
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where V is the Levi-Civita connection of g, and Rg is the corresponding scalar 
curvature, indices are raised and lowered by g. If we specify a Riemannian metric 
5, and a symmetric covariant 2-tensor fc, as initial data on a 3-manifold, they should 
thus in our situation satisfy 

(1.4) Rg - hjk'^ + {iYgkf = 2fM 
and 

(1.5) V.trgfc - = 0, 

because of (|l.3D, ( |l.lD and the fact that u is perpendicular to M. In other words, 
we should also specify the initial value of as part of the data. 

We consider only a restricted class of manifolds M and initial data. The 3-manifold 
M is assumed to be a special type of Lie group, and g, k and /i are assumed to 
be left invariant. In order to be more precise concerning the type of Lie groups 
M = G we consider, let ei, i = 1,2,3 be a basis of the Lie algebra with structure 
constants determined by [ei,ej] = JijCk- If 7j\. = 0, then the Lie algebra and Lie 
group are said to be of class A, and 

(1-6) 7^- = e^.mn'"^ 

where the symmetric matrix n*^ is given by 



(1-7) n^' = 



Definition 1.1. Orthogonal perfect fluid data of class A for Einstein's equations 
consist of the following. A Lie group G of class A, a left invariant Riemannian 
metric g on G, a left invariant symmetric covariant 2-tensor k on G, and a constant 



^0 > satisfying (|1.4[) and (1.5) with ^ replaced by 



We can choose a left invariant orthonormal basis {ci} with respect to g, so that the 
corresponding matrix n*-' defined in (^]^) is diagonal with diagonal elements ni , n2 
and 713. By an appropriate choice of orthonormal basis, ni, ?i2, 713 can be assumed 
to belong to one and only one of the types given in Table |l|. We assign a Bianchi 
type to the initial data accordingly. This division constitutes a classification of the 



class A Lie algebras. We refer to Lemma 21.1 for a proof of these statements 



Let kij — k{ei,ej). Then the matrices rf^ and fc^ commute according to ( |l.5| ), so 
that we may assume kij to be diagonal with diagonal elements ki, k2 and fcs, cf. 



Definition 1.2. Orthogonal perfect fluid data of class A satisfying k2 = k^ and 
^2 — ns or one of the permuted conditions are said to be of Taub type. Data with 
/iQ = are called vacuum data. 

Observe that the Taub condition is independent of the choice of orthonormal basis 



diago naliz in g n a nd k, cf. ( 21.15 ). Considering the equations of Ellis and MacCal- 
lum ( 21.4 )-( |21.8[) , one can see that if n2 = and k2 = k^ at one point in time, 
then the equalities always hold, cf. the construction of the spacetime carried out 
in the appendix. According to , vacuum solutions satisfying these conditions are 
the Taub-NUT solutions. This justifies the following definition. 

Definition 1.3. Taub-NUT initial data are type IX Taub vacuum initial data. 
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Table 1. Bianchi class A. 



Type 


ni 


n2 




I 











II 


+ 








VIo 





+ 




VIIo 





+ 


+ 


VIII 




+ 


+ 


IX 


+ 


+ 


+ 



Definition 1.4. By an orthogonal perfect fluid development of orthogonal perfect 
fluid data of class A, we will mean the following. A connected 4-dimensional Lorentz 



manifold {M,g) and a 2-tensor T, as in (1.1), on {M,g), such that there is an 
embedding i : G M with i*{g) = g, i*{k) — k and i*{p.) — fiQ, where k is the 
second fundamental form of i{G) in {M,g). 

In the appendix, we construct globally hyperbolic orthogonal perfect fluid devel- 
opments, given initial data, and we refer to them as class A developments, cf. 
Definition ^1.1 . We also assign a type to such a development according to the type 



of the initial data. Let us make a division of the initial data according to their 
global behaviour. 

Tiieorem 1.1. Consider a class A development with 1 < 7 < 2. 

1. // the initial data are not of type IX, but satisfy tVgk — 0, then fiQ = and the 
development is causally geodesically complete. Only types I and VIIq permit 
this possibility. 

2. // the initial data are of type I, II, VIq, VIIq or VIII, and satisfy tVgk < 0, 
then the development is future causally geodesically complete and past causally 
geodesically incomplete. Such initial data we will refer to as expanding. 

3. Bianchi IX initial data yield developments that are past and future causally 
geodesically incomplete. Such data are called recollapsing. 

A proof is to be found in the appendix, but observe that this theorem is not new. 
As far as class A developments are concerned, we will restrict our attention to 
equations of state with 1 < 7 < 2. The reason is that there is cause to doubt the 
well posedness of the initial value problem for 2/3 < 7 < 1, cf. |^ p. 85 and p. 88. 
Furthermore, in the Bianchi IX case we use results from p4| concerning recoUapse, 



see Lemma 21.6. In order to be allowed to do that, we need the above mentioned 



condition on 7. What is meant by inextendibility is explained in the following. 

Definition 1.5. Consider a connected Lorentz manifold {M,g). If there is a con- 
nected Lorentz manifold (M, g) of the same dimension, and a map i : M ^ M, 
with i{M) ^ M, which is an isometry onto its image, then {M,g) is said to be 
C^-extendible and {M,g) is called a C^-extension of {M,g). A Lorentz manifold 
which is not C^-extendible is said to be C^-inextendible. 



Remark. There is an analogous definition of smooth extensions. Unless otherwise 
mentioned, manifolds are assumed to be smooth, and maps between manifolds are 
assumed to be as regular as possible. 
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We will use the Kretschmann scalar, 

(1.8) K = Rc.p^sR''^^^ 

as our main measure of whether curvature blows up or not, but in the non- vacuum 
case it is natural to consider the Ricci tensor contracted with itself Rc/^R"^ . The 
next theorem states the main conclusion concerning developments. 

Theorem 1.2. For class A developments with 1 < 7 < 2, we have the following 
division. 

1. Consider expanding initial data of type I, II or VIIq with 1 < 7 < 2 which are 
not of Taub vacuum type. Then the Kretschmann scalar is unbounded along 
all inextendible causal geodesies in the incomplete direction. 

2. Consider non-Taub-NUT recollapsing initial data with 1 < 7 < 2. Then the 
Kretschmann scalar is unbounded along all inextendible causal geodesies in 
both incomplete directions. 

3. Expanding and recollapsing data with 7 = 2 and fiQ > 0. Then the Kretschmann 
scalar is unbounded along all inextendible causal geodesies in all incomplete 
directions. 

4. Expanding and recollapsing data with /.iq > 0. Then RapR"^ is unbounded 
along all inextendible causal geodesies in all incomplete directions. 

In all cases mentioned above the class A development is -inextendible. 

Remark. Observe that the Bianchi VIII vacuum case was handled in and the 
Bianchi VIq vacuum case in [^] . The above theorem thus isolates the vacuum Taub 
type solutions as the only ones among the Bianchi class A spacetimes that do not 
exhibit curvature blow up, given our particular matter model. 

We now turn to the results that are expressed in terms of the variables of Wainwright 
and Hsu. The equations and some of their properties are to be found in Section |[ 
The appendix contains a derivation. It is natural to divide the matter models into 
two categories; the non-stiff fluid case and the stiff fluid case (7 = 2). 

Let us begin with the non-stiff fluid case, including the vacuum case. We confine 
our attention to Bianchi IX solutions. The existence interval stretches back to 
—00 which corresponds to the singularity. There are some fixed points to which 
certain solutions converge, and data which lead to such solutions together with 
data of Taub type will be considered to be non-generic. The Kasner map, which 
is supposed to be an approximation of the Bianchi IX dynamics as one approaches 
a singularity, is illustrated in Figure |l} The circle in the S_|_E_-plane appearing 
in the figure is called the Kasner circle, and we have depicted two bounces of the 
Kasner map. The starting point is marked by a star, and the end point by a plus 
sign. Given a point x on the Kasner circle, the Kasner map yields a new point 
y on the Kasner circle by taking the corner of the triangle closest to x, drawing 
a straight line from the corner through x, and then letting y be the second point 
of intersection between the line and the Kasner circle. One solid line corresponds 
to the closure of a vacuum type II orbit of the equations of Wainwright and Hsu. 
Actually, it is the projection of the closure of such an orbit to the I]+S_-plane. 
A vacuum type II solution has one Ni non-zero and the other zero, and the three 
different Ni correspond to the three corners of the triangle; the rightmost corner 
corresponds to iVi 7^ and the corner on the top left corresponds to N3 ^ 0. The 
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-1.5 -1 -0.5 



0.5 1 1.5 



Figure 1. The Kasner map. 



constraint ( |2.3| ) for the vacuum type II solutions is given by 

^9 . ^9 3 



1. 



The closure of this set is given a name in the following definition. 
Definition 1.6. The set 

A = {{n, S+, iVi, iV2, iVs) : + \NiN2\ + lA^aA^al + = 0} n M, 



where M is defined by (2.3), is called the Bianchi attractor. 



The main result of this paper is that for generic Bianchi IX data, the solution 
converges to the attractor. That is 

(1.9) lim (f7 + iViiV2 + 7V2A^3 + iV3iVi) = 0. 



r — > — oo 



This conclusion supports the statement that the Kasner map approximates the 
dynamics, and also the statement that the matter content loses significance close 
to the singularity. Let us introduce some terminology. 

Definition 1.7. Let / e C°°(M",M"), and consider a solution x to the equation 

^ = fox, x{0) = xo, 

with maximal existence interval (i_, t-|-). We call a point an a-limit point of the 
solution X, if there is a sequence tk t- with x{tk) x^. The a-limit set of x 
is the set of its a-limit points. The w-limit set is defined similarly by replacing t_ 
with 

Remark. If i_ > —oo then the a-limit set is empty, cf. [ p^ . 

Thus, the a-limit set of a generic solution is contained in the attractor. The desired 
statement is that the a-limit set coincides with the attractor, but the best result we 
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have achieved in this direction is that there must at least be three a-hmit points on 
the Kasner circle. This worst case situation corresponds to the solution converging 
to a periodic orbit of the Kasner map with period three. Observe that we have not 
proven anything concerning Bianchi VIII solutions. 

Let us sketch the proof. It is natural to divide it into two parts. The first part 
consists of proving the existence of an a-limit point on the Kasner circle. We 
achieve this in the following steps. First we analyze the a-limit sets of the Bianchi 
types I, II and VIIq. An analysis of types I of II can also be found in Ellis and 
Wainwright jl^ . Then we prove the existence of an a-limit point for a generic 
Bianchi IX solution. To go from the existence of an a-limit point to an a-limit 
point on the Kasner circle, we use the analysis of the lower Bianchi types. In the 
second part, we prove ( |1.9| ). Let d be the function appearing in that equation. We 
assume that d does not converge to zero in order to reach a contradiction. The 
existence of an a-limit point on the Kasner circle proves that there is a sequence 
Tk — oo such that d(rfc) ^0. If d does not converge to zero there is a 5 > 0, and a 
sequence Sk —oo such that d{sk) > S. We can assume Sk < Tk and conclude that 
d on the whole has to grow (going backwards) in the interval [sk,Tk]. What can be 
said about this growth? In Section ^ we prove that we can control the density 
parameter in this process, assuming S is small enough, which is not a restriction. 
As a consequence can be assumed to be arbitrarily small during the growth. 



Some further arguments, given in Section 15, show that we can assume the growth 



to occur in the product N2N3, using the symmetries of the equations. Furthermore, 
one can assume the -variables to be arbitrarily close to = (—1,0), 

and that some expressions dominate others. For instance 1-I-E+ can be assumed to 
be arbitrarily much smaller than N2N3. This control introduces a natural concept 
of order of magnitude. The behaviour of the product N2N3 will be oscillatory; it 
will look roughly like a sine wave. The point is to prove that the product decays 
during a period of its oscillation; that would lead to a contradiction. The variation 
during a period can be expressed in terms of an integral, and we use the order of 
magnitude concept to prove an estimate showing that this integral has the right 
sign. 

Now consider the stiff fluid case with positive density parameter. In this case we will 
consider Bianchi VIII and IX solutions. The analysis is similar for the other cases 
and a description of the results is to be found in Section |l^. Again the singularity 
corresponds to —00. The density parameter converges to a non-zero value, all 
the Ni converge to zero, and in the E-(_I]_-plane the solution converges to a point 
inside the triangle shown in Figure ^. 

In Section H, we formulate the equations of Wainwright and Hsu and briefly describe 
their origin and some of their properties. Section || contains some elementary 
properties of solutions. We give the existence intervals of solutions to the equations, 
and prove that the fiS+S- -variables are contained in a compact set to the past 
for Bianchi IX solutions. As in the vacuum case, we also prove that (S-|_, S_) can 
converge to (—1, 0) only if the solution is of Taub type, although this is no longer 
a characterization. In Section ^, we mention some critical points and make more 
precise the statement that solutions converging to these points are non-generic. 
Included in this section are also two technical lemmas relevant to the analysis. The 
monotonicity principle is explained in Section ^. It is fundamental to the analysis 
of the a-limit sets of the solutions. We present two applications; the fact that all 
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-1.5 -1 -0.5 



0.5 1 1.5 



Figure 2. The triangle mentioned in the text. 



a-hmit points of Bianchi IX solutions are of type I, II or VIIo and an analysis of the 
vacuum type II orbits. The last application is not complicated, but illustrates the 
arguments involved as well as demonstrating how the map depicted in Figure [l] can 
be viewed as a sequence of type II orbits. Section |^ deals with situations such that 
one has control over the shear variables and the density parameter. Specifically, 
it gives a geometric interpretation of some of the equations in riS-|-I]_-space. As 
an application, we prove that if a Bianchi IX solution has an a-limit point on the 
Kasner circle then all the points obtained by applying the Kasner map to this 
point belong to the a-limit set of the solution. The stiff fluid case is handled 
in Section!^ In this case the a-limit set consists of a point regardless of type. 



Sections pi- 10 deal with the lower order Bianchi types needed in order to analyze 



Bianchi IX. An analysis of types I of II can also be found in Ellis and Wainwright 
18). Section |ri| gives the possibilities for a Taub type Bianchi IX solution. The 
technical Section 12 is needed in order to prove the existence of an a-limit point 



for Bianchi IX solutions, and also to prove that the set of vacuum type II points 
is an attractor. It is used for approximating the solution in situations where the 
behaviour is oscillatory. Section |l^ proves the existence of an a-limit point for a 
Bianchi IX solution and the existence of an a-limit point on the Kasner circle for 
generic Bianchi IX solutions. In Section ^ we prove that if one has control over 
the sum |iViA^2| + 1-^2-^31 + lA'^sA^i] in some time interval [Ti,r2], and control over 
O in T2 then one has control over SI in the entire interval. This rather technical 
observation is essential in the proof that generic solutions converge to the attractor. 
The heart of this paper is Section |l^ which contains a proof of ( 1 . 9 ) . It also contains 
arguments that will be used in Section ^ to analyze the regularity of the set of 
non-generic points. In Section |l^, we observe that the convergence to the attractor 
is uniform, and in Section |l^ we prove the existence of at least three non-special 
a-limit points on the Kasner circle. We formulate the main conclusions and prove 
Theorem 1.2 in Sectio n [l9| . In Section |2^, we relate our results concerning stiff fluid 
solutions to those of 0. The appendices contain results relating solutions to the 
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equations of Wainwright and Hsu with properties of the class A developments and 
some curvature computations. 



2. Equations of Wainwright and Hsu 
The essence of this paper is an analysis of the asymptotic behaviour of solutions 



to the equations of Wainwright and Hsu (2.1)-(2.3). One important property of 
these equations is that they describe all the Bianchi class A types at the same time. 
Another important property is that it seems that the variables remain in a compact 
set as one approaches a singularity. In the Bianchi IX case, this follows from the 
analysis presented in this paper. Let us give a rough description of the origin of the 
variables. In the situations we consider, there is a foliation of the Lorentz manifold 
by homogeneous spacelike hypersurfaces diffeomorphic to a Lie group G of class A. 
One can define an orthonormal basis Bq, a = 0, 3, such that ej, i = 1,2, 3, span 
the tangent space of the spacelike hypersurfaces of homogeneity, and eo = dt for 
a suitable globally defined time coordinate t. It is possible to associate a matrix 



with the spacelike vectors e^, as in (1.7), and assume it to be diagonal with 



diagonal components n^. One changes the time coordinate by dt/dr = 3/0, where 
is minus the trace of the second fundamental form of the spacelike hypersurface 
corresponding to t. The Ni{T) below are the ^^(r) divided by 6{t), the E+ and I]_ 
correspond to the traceless part of the second fundamental form of the spacelike 
hypersurface corresponding to r, similarly normalized, and finally Vl = 3^/9^. We 
will refer to E+ and S_ as the shear variables, and to f2 as the density parameter. 
The question then arises to what extent this makes sense, since 9 could become 
zero. An answer is given in the appendix. For all the Bianchi types except IX, this 
procedure is essentially harmless, and the variables of Wainwright and Hsu capture 
the entire Lorentz manifold. In the Bianchi IX case, there is however a point at 
which = 0, at least if 1 < 7 < 2, see the appendix, and the variables are only 
valid for half a development in that case. As far as the analysis of the asymptotics 
are concerned, this is however not important. A derivation of the equations is given 
in the appendix. They are 



(2.1) 





- (9- 


4S+)7Vi 






= (9 + 


2E+ + 2V3I]_ 






= (9 + 


2E+ - 2V3I]- 


)N3 




= -(2 


- q)^+ - 35+ 






= -(2 


- (7)E_ - 35- 




Vt' 


= [2<Z- 


-(37-2)]r!. 





The prime denotes derivative with respect to a time coordinate r, and 

q = i(37-2)r! + 2(S^ + S^) 
(2.2) 5+ = ^[{N2 - N^f - Nii2Ni - N2 - N3)] 



10 



HANS RINGSTROM 



The constraint is 

(2.3) n + j:l + j:l + ^[n^ + iv| + nI - 2{NiN2 + N2N3 + N3N1)] ^ 1. 



We demand that 2/3 < 7 < 2 and 17 > 0. The equations (2.1)-(2.3) have certain 
symmetries, described in Wainwright and Hsu jl^. By permuting Ni, N2, N3 arbi- 
trarily, we get new solutions, if we at the same time carry out appropriate combina- 
tions of rotations by integer multiples of 27r/3, and reflections in the (E_(-, I]_)-plane. 
Explicitly, the transformations 

(7Vi,iV2,7V3) = {N3,N,,N2), (S+,E_) = (-is+ + ^a/SE-, -^^32+ - ^S-) 
and 

(7Vi,iV2,7V3) = (iVi,7V3,iV2), (S+,S_) - (S+,-S-) 

yield new solutions. Below, we refer to rotations by integer multiples of 27r/3 
as rotations. Changing the sign of all the Ni at the same time does not change 
the equations. Classify points Ni, N2, N3) according to the values of 

Ni,N2,N3 in the same way as in Table |l|. Since the sets Ni > 0, Ni < and 
Ni — are invariant under the flow of the equations, we may classify solutions to 



(Ej)-(2-3) accordingly. 



Definition 2.1. The Kasner circle is defined by the conditions Ni = fl = and 
the constraint ( |2.3| ). There are three points on this circle called special: (S+, E_) = 
(-1,0) and (l/2,±V3/2). 

The following reformulation of S'^ is written down for future reference, 
(2.4) 

= -(2 - 20 - 2E^ - 2S2_)(I]+ + 1) - ^(2 - 7)0^+ + ^iVi(iVi - N2 - N3). 



3. Elementary properties of solutions 

Here we collect some miscellaneous observations that will be of importance. Most of 
them are similar to results obtained in . The a-limit set defined in Definition 
plays an important role in this paper, and here we mention some of its properties. 

Lemma 3.1. Let f and x he as in Definition The a-limit set of x is closed 
and invariant under the flow of f . If there is a T such that x{t) is contained in a 
compact set for t < T, then the a-limit set of x is connected. 

Proof See e. g. □ 

Definition 3.1. A solution to ( |2.lD -( ^^ satisfying N2 = N^ and S_ = 0, or one 

of the conditions found by applying the symmetries, is said to be of Taub type. 

Remark. The set defined by N2 = N3 and E_ = is invariant under the flow of 



(2.1) 



Lemma 3.2. The existence intervals for all solutions to (2.1 )-{2.c ) except Bianchi 



IX are (—00,00). For Bianchi IX solutions we have past global existence. 
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Proof. As in the vacuum case, see IQ. □ 

By observations made in the appendix, — oo corresponds to the singularity. 

Lemma 3.3. Let 2/3 < 7 < 2. Consider a solution of type IX. The image 
00, 0]) is contained in a compact set whose size depends on the initial 
data. Further, if at a point in time N3 > N2 > Ni and N3 > 2, then N2 > N^/IO. 

Proof. As in the vacuum case, see □ 

That S_, fl) is contained in a compact set for aU the other types foUows from 
the constraint. The second part of this lemma will be important in the proof of 
the existence of an a-limit point. One consequence is that one Ni may not become 
unbounded alone. 

The final observation is relevant in proving curvature blow up. One can define a 



normalized version (22.3) of the Kretschmann scalar (1.8), and it can be expressed 
as a polynomial in the variables of Wainwright and Hsu. One way of proving that a 
specific solution exhibits curvature blow up is to prove that it has an a-limit point 
at which the normalized Kretschmann scalar is non-zero. We refer to the appendix 
for the details. It turns out that this polynomial is zero when N2 = N3, Ni — 0, 
I]_ = 0, S+ = — 1 and = 0. The same is true of the points obtained by applying 
the symmetries. It is then natural to ask the question: for which solutions does 
(S+,S_) converge to (-1,0)? 



Proposition 3.1. A solution to (2.1 )- ^2.3l ) with 2/3 < 7 < 2 satisfies 

lim (E+(T),S_(r)) = (-1,0), 

r — * — 00 

only if it is contained in the invariant set I]_ = and N2 = N3. 



Remark. The proposition does not apply to the stiff fluid case. The analogous 
statements for the points (S+,I]_) = (1/2, ±-\/3/2) are true by an application 
of the symmetries. We may not replace the implication with an equivalence, cf. 



Proposition 9.1 



Proof. The argument is essentially the same as in the vacuum case, see |l6|. We 
only need to observe that Q will decay exponentially when (I]_|-,S_) is close to 
(-1,0). □ 



4. Critical points 



Definition 4.1. The critical point F is defined hy fl — 1 and all other variables 
zero. In the case 2/3 < 7 < 2, we define the critical point P^{II) to be the type II 
point with S_ = 0, A^i > 0, S+ = (37 - 2)/8 and = 1 - (37 - 2)/16. The critical 
points P^{II), i = 2,3 are found by applying the symmetries. 

It will turn out that there are solutions which converge to these points as r ^ — 00. 
The main objective of this section is to prove that the set of such solutions is small. 
Observe that only non-vacuum solutions can converge these critical points. 

Definition 4.2. Let lyiio denote initial data to ( |2.l| )-( |273| ) of type VIIq with SI > 0, 
and correspondingly for the other types. Let T^viio be the elements of lyiio such that 
the corresponding solutions converge to one of P^{II) as r ^ — 00 and similarly 
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for Bianchi II and IX. Finally, let Twig be the elements of Iviio such that the 
corresponding solutions converge to as r — s- — oo, and similarly for the other 
types. 

Remark. The sets T\i and so on depend on 7, but we omit this reference. 

Observe that Xj, Xn, IviIq and Xix are submanifolds of of dimensions 2, 3, 4 
and 5 respectively. They are diffcomorphic with open sets in a suitable M"; project 
Q, to zero. We will prove that V\\ consists of points and that Ti is the point F. 



Let 2/3 < 7 < 2 be fixed. In Theorem 16.1, we will be able to prove that the sets 
•^iIj-^ViIqj -^iXj T-'yUo and T^ix are submanifolds of of dimensions 1, 2, 3, 1 
and 2 respectively. This justifies the following definition. 

Definition 4.3. Let 2/3 < 7 < 2. A solution to (p.l|)-(p7^) is said to be generic if 
it is not of Taub type, and if it does not belong to Ti, Tn, Tyu^ , .Fix, ^11, T^viIq or 

We will need the following two lemmas in the sequel. 



Lemma 4.1. Consider a solution x to (2.1)-(2.S) such that x has P^{II) as an 



a-limit point but does not converge to it. Then x has an a-limit point of type II, 
which is not P^{II). 

Remark. There is no solution satisfying the conditions of this lemma, but we will 
need it to establish that fact. 

Proof. Consider the solution to belong to R^, and let the point represent P^{II). 
There is an e > such that for each T, there is a r < T such that x{t) does not 
belong to the open ball B^{xq). In one can compute that 

q + 2S+ ± 2V3I]_ > 0. 

Let e be so small that these expressions are positive in B^{xq). Let r^, — > — cx) 
be a sequence such that x(Tk) Xq, and let < tj^ be a sequence such that 
x{sk) € OB^^xq) and x{{sk,Tk]) C B^{xo). Since x{sk) is contained in a compact 
set, there is a convergent subsequence yielding an a-limit point which is not P^ {II). 
Since N2 and converge to zero in r^. and decay in absolute value from r^, to s^,, 
the a-limit point has to be of type II {Ni has to be non-zero for the new a-limit 
point if e is small enough). □ 



Lemma 4.2. Consider a solution x to (2.1)-(2.o) such that x has F as an a-limit 



point, but which does not converge to F. Then x has an a-limit point of type I 
which is not F. 



Remark. The same remark as that made in connection with Lemma LI holds 
concerning this lemma. 

Proof. The idea is the same as the previous lemma. We need only observe that 
q - 41]+, q + 2E+ + 2^/31]- and q + 2I]+ - 2\/3E- are positive in F. □ 

5. The monotonicity principle 



The following lemma will be a basic tool in the analysis of the asymptotics, we will 
refer to it as the monotonicity principle. 
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Lemma 5.1. 

(5.1) 



Consider 



dx 
'dt 



fox 



where f G C°°(M",IE 
invariant under the 



"'). Let U be an open subset of I 
low of the vectorfield f. Let G 



" , and M a closed subset 
U ^ R &e a continuous 



function such that G{x{t)) is strictly mon oton e for any solution x{t) of (SA), as 
long as x{t) € UHM . Then no solution of (5.1 ) whose image is contained in U flM 
has an a- or u-limit point in U . 

Remark. Observe that one can use M — R " . W e will mainly choose AI to be 
the closed invariant subset of R^ defined by ( |2.3| ). If one Ni is zero and two are 
non-zero, we consider the number of variables to be four etc. 

Proof. Suppose p e L/ is an a- limit point of a solution x contained in UO M . Then 
Go a; is strictly monotone. There is a sequence tn — > t^ such that x{tn) — > p by our 
supposition. Thus G{x{tn)) G{p), but Gox is monotone so that G{x{t)) G{p). 
Thus G{q) — G{ p) for all a-limit points q of x. Since M is closed p G M . The 
solution a; of (5.1), with initial value p, is contained in M by the invariance property 
of M, and it consists of a-limit points of x so that G{x{t)) = G{p) which is constant. 
Furthermore, on an open set containing zero it takes values in U contradicting the 
assumptions of the lemma. □ 

Let us give an example of an application. 

Lemma 5.2. Consider a solution to (2.1)-(2.j) of type VIII or IX. If it has an 

a-limit point, then 

{NiN2N3){t) =0. 



lim 

" — ^ — oo 



Proof. Let U of Lemma 5.1 be defined by the union of the sets Ni ^ 0, i 
M by the constraint (2.3), and G by the function NiN2N3. Compute 

(5.2) (TViiVaTVs)' = 3qNiN2N3. 



1,2,3, 



Consider a solution x of (2.1)-(2.3). We need to prove that Gox is strictly monotone 
as long as x{t) e C/ n M. By (^.2|) the only problem that could occur is q = 0. 
However, q = implies jS'^j -I- > by (2.1)-(2.3) so that Go a; has the desired 
property. If the sequence Tk —oo yields the a-limit point we assume exists, then 
we conclude that 

Since NiN2Nr^ is monotone, we conclude that it converges to zero. □ 

One important consequence of this observation is the fact that all a-limit points 
of Bianchi VIII and IX solutions are of one of the lower Bianchi types. Since the 
a-limit set is invariant under the flow, it is thus of interest to know something about 
the a-limit sets of the lower Bianchi types, if one wants to prove the existence of 
an a-limit point on the Kasner circle. 

Let us now analyze the vacuum type II orbits and define the Kasner map. 



Proposition 5.1. A Bianchi II vacuum solution of (2.1)-(2.c) with Ni > and 
N2 — N3 = satisfies 



(5.3) 



lim iVi = 0. 

"— >±oo 
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The oj-limit set is a point in /Ci and the a-limit set is a point on the Kasner circle, 
in the complement of the closure of K-i . 



Remark. What is meant by Ki is explained in Definition 6.1 



Proof. Using the constraint (2.3) we deduce that 

We wish to apply the monotonicity principle. There are three variables. Let U be 



defined by TVi > 0, M be defined by (|2J), and G(S+, S_, iVi) = E+. We conclude 
tha t (^ ) is true as follows. Let r„ — ^ cx3. A subsequence yields an w-limit point 
by (|2.3|). The monotonicity principle yields iVi(r„j^) for the subsequence. The 



argument for the a-limit set is similar, and equation (5.3) follows. Combining this 
with the constraint, we deduce 

lim q = 2. 

Using the monotonicity of S+, we conclude that (S+,I]_) has to converge. As 
for the a-limit set, convergence to /Ci is not allowed since < close to /Ci. 
Convergenc e to one of the special points in the closure of /Ci is also forbidden, since 
Proposition |3.l| would imply A^i = for the solution in that case. Assume now that 
(S+,I]_) ((T+,cr_) as r ^ CO. Compute 



We get 

for arbitrary (E+,S_) belonging to the solution. Since = (q — 4E+)A^i and 
iVi 0, we have to have cr+ > 1/2. If a+ = 1/2, then a- = ±V3/2. The 
t wo corresponding lines in the I]_|_I]_-plane, obtained by substituting (tT+, ct-) into 
(pTsI), do not intersect any points interior to the Kasner circle. Therefore ct+ — 1/2 
is not an allowed limit point, and the proposition follows. □ 

Observe that by ( ^.4|) , the projection of the solution to the S+E_-plane is a straight 
line. The orbits when A''2 > and when iVa > are obtained by applying the 
symmetries. Figure shows a sequence of vacuum type II orbits projected to the 
S+E_-plane. The first line, starting at the star, has A^i > 0, the second > 
and the third N2 > 0. 

Definition 5.1. If xq is a non-special point on the Kasner circle, then the Kasner 
map applied to xq is defined to be the point xi on the Kasner circle, with the 
property that there is a vacuum type II orbit with xq as an tj-limit point and xi 
as an a-limit point. 



6. Dependence on the shear variables 



In several arguments, we will have control over the shear variables and the density 
parameter in some time interval, and it is of interest to know how the remaining 
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variables behave in such situations. Consider for instance the expression muhiplying 
A^i in the formula for N{, see (^]^). It is given by g — 4I]-|_ and equals zero when 

(6.1) i(37-2)f]+(l-S+)2 + I]2_ ^1. 

The set of points in f7S_|_S_-space satisfying this equation is a paraboloid, and 
the intersection with = is the dashed ci rcle shown in Figure ^. If {fl, S+, E_) 
belongs to the interior of the paraboloid ( |6.l[ ) with f2 > 0, then |iVi|' will be 
negative, so that \Ni \ increases as we go backward. Outside of the paraboloid, |A^i| 
decreases. The situation is similar for N2 and N3. Observe that the circle obtained 
by letting 57 = in ( |6.l[ ) intersects the Kasner circle in two special points. The 
same is true of the rotated circles corresponding to N2 and N3. It will be convenient 
to introduce notation for the points on the Kasner circle at which \Ni\' is negative. 



1.5 - 



1 - 
0.5 - 

w' - 




-1.5 - 



-2' ' ' ' ' ' ' ' ' 

-1.5 -1 -0.5 0.5 1 1.5 2 2.5 



Figure 3. The circles mentioned in the text. 

Definition 6.1. Wc let ICi, IC2 and /C3 be the subsets of the Kasner circle where 
q - 41]+ <0, q + + 2x/3i;_ < and g + 2S+ - 2\/3E_ < respectively. 

Remark. On the Kasner circle, = so that q = 2(I]+ -|- S^) = 2 under the 
conditions of this definition. 

It also of interest to know when the derivatives of N2N3 and similar products are 
zero. Since (iV2A^3)' = (2g -I- 4E+)7V2-^3, we consider the set on which q + 2E+ 
equals zero. This set is a paraboloid and is given by 

1(37 -2)f} + (S+ + + = 

The intersection with the plane f2 = is the circle with radius 1 /2 shown in Figure 
|[ Again, inside the paraboloid |A^2-^3| increases as we go backward, and outside it 
decreases. There are corresponding paraboloids for the products N1N2 and N1N3. 
Observe that in the non-vacuum case, it is harmless to introduce lj = 17^/^ and 
then the paraboloids become half spheres. 
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Proposition 6.1. Consider a Bianchi IX solution to (2.1)- (2. 5) with 2/3 < 7 < 2. 
// the solution has a non-special a-limit point x on the Kasner circle, then the 
closure of the vacuum type II orbit with x as an Lu-limit point belongs to the a-limit 
set. 

Remark. The same conclusion holds for a Bianchi type VIIq solution with A^i — 0, 
if it has an a-limit point in IC2 or /C3. 

Proof. Assume the limit point lies in /Ci with = ((t+,(T_). There is a 

sequence — > —00, such that the solution evaluated at converges to the point 
on the Kasner circle. There is a ball _B^(CT+,cr_) in the I]+I]_-plane, centered at 
this point, such that |A^2|, l-^ali |-^i-^2|, l-^i-^al and fl all decay exponentially, at 
least as e^'^ for some fixed ^ > 0, and A^i increases exponentially, at least as e^^'^ , in 
the closure of this ball. There is a K such that (E+(Tfc), E_(Tfe)) G Bji{a+, (t_) for 
all k > K. For each time we enter the ball, we must leave it, since if we stay in it to 
the past, A^i will grow to infinity whereas N2 and N3 will decay to zero, in violation 
of the constraint. Thus for each Tk, k > K, there is a ifc < corresponding to the 
first time we leave the ball, starting at and going backward. We may compute 



where 

\h{T)\ < efce^f^"^") 

in [tk, Tk] and ^ 0. Thus 

E_(Tfc) S_fa) 
2-S+(Tfe) 2-S+(ffe) 



hdr. 



But 



and in consequence 



hdr] < 



E_(rfc) S_(4 



0. 



2-I]+(Tfc) 2-^+itk) 
We thus get a type II vacuum limit point with A^i > 0, to which we may apply the 
flow, and deduce the conclusion of the lemma. The statement made in the remark 
follows in the same way. Observe that the only important thing was that the limit 
point was in /Ci and Ni was non-zero for the solution. □ 

7. The stiff fluid case 

In this section we will assume $7 > and 7 = 2 for all solutions we consider. We 
begin by explaining the origin of the triangle shown in Figure |[ Then we analyze 
the type II orbits. They yield an analogue of the Kasner map, connecting two 



points inside the Kasner circle, and we state an analogue of Proposition 3.1 for this 
map. We then prove that is bounded away from zero to the past. Only in the 
case of Bianchi IX is an argument required, but this result is the central part of 
the analysis of the stiff fluid case. A peculiarity of the equations then yields the 
conclusion that \NiN2\ -\- |A^2A^3| + l-^a-^il converges to zero exponentially. This 
proves that any solution is contained in a compact set to the past, and that all 
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Qf-limit points are of type I or II. Another consequence is that has to converge to 
a non-zero value; this requires a proof in the Bianchi IX case. Next one concludes 
that all Ni converge to zero, since if that were not the case, there would be an 
a-limit point of type II to which one could apply the flow, obtaining a-limit points 
with different 17:s. Then if a Bianchi IX solution had an a-limit point outside the 
triangle, one could apply the 'Kasner' map to such a point, obtaining an a-limit 
point with some iVj > 0. Finally, some technical arguments finish the analysis. 

In the case of a stiff fluid, that is 7 = 2, it is convenient to introduce 

We then have, since 87 — 2 = 4, 

(7.1) uj' ^ -(2-q)uj. 

The expression ^l+E^ + E^ turns into lj^ + S'^ + S^, and the uj, S+, E_ -coordinates 
of the type I points obey 

(7.2) Lu^ + T,l + Y.l^l, uj>0. 

In the stiff fluid case, all the type I points are fixed points, and they play a role 
similar to that of the Kasner circle in the vacuum case. 

Let us make some observations. If iVi ^ 0, then A^^ = is equivalent to 5— 4E+ = 0. 
Dividing by 2 and completing squares, we see that this condition is equivalent to 

(7.3) uj"^ + + -£1=^1, uj>0. 

By applying the symmetries, the conditions — 0, Ni ^ are consequently 
all fulfilled precisely on half spheres of radii 1. Since \Ni\' < corresponds to 
an increase in |iVi| a s we go backward, |A^i| increases exponentially as we are 



inside the half sphere ( 7.3) an d decreases exponentially as we are outside it. I f on e 



takes the intersection of ( fr.2|) and ([7.3|), one gets the subset S_|_ = 1/2 of ( |7.2D 

The corresponding intersections for A^2 and A3 yield two more lines in the 

plane. T oget her they yield the triangle in Figure Consequently, if (w, E+, E_) is 



close to (7.2) and (S+, S_) is in the interior of the triangle, then all the Ni decay 
exponentially as t — ^ —00. 

Let Ml be the subset wS+S.-space obeying ( [7.2|) with a; > and E+ > 1/2 and 
M2, M3 be the corresponding sets for N2 and N3. We also let Ci be the subset 



of the intersection between ( |7.2[) and ( |7.3[) with a; > and correspondingly A^2 and 
A^3 yield £2 and £3. 



Lemma 7.1. Consider a solution to (2.1)- ^2.q ) with 7 = 2 such that Ni > 0, 
LU > and N2 = N3 = 0. Then 

(7.4) lim Ni{t) = 

r — *±ckj 



and (w, E+, E_) converges to a point, satisfying (7J ) anduj > 0, in the complement 



of Ci U Ml, as T ^ —00. In ujT,+Ti -- sp ace, the orbit of the solution is a straight 



line connecting two points satisfying (TJj. Ifuj>0, it is strictly increasing along 
the solution, going backwards in time. 

Proof. Since q < 2 for the entire solution, we can apply the monotonicity principle 



with U defined by g < 2, G defined by and M by the constraint (2.3). If q 



does not converge to 2 as r —> —00, we get an a-limit point with q < 2. We have 
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a contradiction. This argument also yields the conclusion that A^i 
Equation (7.4) follows. Observe that 



as r ^ oo. 



(7.5) 

and 

(7.6) 



Consequently, 
T — !■ CX3 and as r 



S]_ and are all monotone so that they converge, both as 
- —oo. It also follows from (7.5) and ( [7.6[ ) that the quotients 
(2 — I]+)/a; and are constant. Thus the orbit in cjS]+S„-space describes 

a straight line connecting two points satisfying (7.2). As r ^ — oo, the solution 
cannot converge to a point in £i U A^i for the following reason. Assume it does. 
Since S+ decreases as r decreases, see (7.5), we must have E_|_ > 1/2 for the entire 
solution, since by assumption converges to a value > 1/2. But then < 
for the entire solution by (2.1) and (2.3). Thus iVi increases as we go backward, 
contradicting the fact that A^i ^ 0. □ 

The next thing we wish to prove is that if a solution has an a-limit point x in the 
set A4i, and A^i 7^ for the solution, then we can apply the 'Kasner' map to that 
point. What wc mean by that is that an entire type II orbit with x as an w-limit 
point belongs to the a-limit set of the original solution. From this one can draw 
quite strong conclusions. Observe for instance that by (7.1), lu is monotone for 
a Bianchi VIII solution to (2.1)-(2.3). Thus uj converges as r —00 since it is 
bounded. If the Bianchi VIII solution has an a-limit point of type I outside the 
triangle, we can apply the Kasner map to it to obtain a-limit points with different 
UJ. But that is impossible. 



Lemma 7.2. Consider a solution to (2.1)-(2.i) with 7 = 2 such that Ni ^ 0. 
Then if the solution has an a-limit point x G Mi, the orbit of a type II solution 
with X as an Lu-limit point belongs to the a-limit set of the solution. 



Proof. The proof is analogous to the proof of Proposition 6.1 



□ 



Consider a solution such that uj > 0. We want to exclude the possibility that uj ^ 
as T ^ —00. Considering (7.1), we see that the only possibility for uj to decrease is 
if q > 2. In that context, the following lemma is relevant. 



Lemma 7.3. Consider a Bianchi IX solution to (2.1)-(2.o) with 7 = 2. There is 
an ao such that if a < ao and 

iNiN2N3){T) < a, 

then 

q{T) - 2 < 4ai/^ 

Proof. By a permutation of the variables, we can assume iVi < N2 < -/V3 in r. 
Observe that 

g-2<37Vi(7V2+iV3) 
by the constraint (2.3). If A^3 < a^/^ in r, we get g — 2 < 6a < 4a^/^ if ao is small 



enough. If > a^^^^in r, we get 



N1N2 < a^/^ 
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Assume, in order to reach a contradiction, {NiN3){t) > a^/^^. Then N2{t) < c?!'^ ^ 
so that iVi(r) < c?!'^ and Nzij) > oT^^^ . By Lemma 3^ we get a contradiction if 
ao is smaU enough. Thus 

g - 2 < ■i{NYN2 + iViiV3)(r) < 3(a^/3 + a^l'^) < 

if ao is smah enough. □ 

For aU solutions except those of Bianchi IX type, u is monotone increasing as r 
decreases. Thus, uj is greater than zero on the a-hmit set of any non-vacuum 
solution which is not of type IX. It turns out that the same is true for a Bianchi 
IX solution. 



Lemma 7.4. Consider a Bianchi IX solution to (2.1)-(2.S) with 7 = 2 such that 
LU > 0. Then there is an e > such that lo^t) > e for all t < 0. 

Proof. Assume all the Ni are positive. The function 

9 = 

satisfies cj)' = 2(j). Thus, for t < 0, 

(iViiV27V3)'/'(r) = a;(T)0(O)e2- < Ce^-, 



because of Lemma 3.3. For r < T < 0, we can thus apply Lemma 7.3, so that for 
T<T, 



iq{s) - 2)ds = J {q{s) - 2)ds + J {q{s) - 2)ds < 4C J e^'ds+ 

/ {q{s) - 2)ds < 2Ce^^ + / (g(s) - 2)ds < C < 
Jt Jt 



Consequently, 



uj{t) = Lu{0) cxp(- / {q{s) - 2)ds) > uj{Q)e 



-C 



and the lemma follows. □ 

The next lemma will be used to prove that uj converges for a Bianchi IX solution. 



Lemma 7.5. Consider a solution to (2.1 )-{2.c ) with 7 = 2 and a; > 0. Then there 
is an a > Q and a T such that 



\NiN2\ + liVziVal + iTVaA^il < e"" 



for all T <T . 



Proof Consider g = iTVaA^al/^^- Then 

g' = (2u;2 + 2(1 + Y.+f + 2Y?_)g. 
Since w(r) > e for all r < 0, we conclude that 

ff(r) < ff(0) exp(2e2T) 

so that 

|(iV2A^3)(T)| <ff(0V(T)exp(2eV). 
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There are similar estimates for the other products. By Lemma 3.3, we know that 
uj is bounded in (—00, 0] so that by choosing a — and T negative enough the 
lemma follows. □ 



Corollary 7.1. Consider a solution to (2.1)-(2.5) with 7 = 2 and uj > 0. Then 
(w, S_, A^i, iV2, A^3)(— 00, 0] is contained in a compact set and all the a-limit 
points are of type I or II. 



Lemma 7.6. Consider a solution to (2.1)- ^.3^ with 7 = 2 and ui > Q. Then 

lim lo{t) = cjo > 0. 



r — > — 00 



Proof. Since this follows from the monotonicity of io in all cases except Bianchi IX, 
see ([7.l[ ), we assume that the solution is of type IX. Let Tfc — > — 00 be a sequence 
such that u!{Tk) cji > 0. This is pos sible since cu is constrained to belong to 
a compact set for r < by Lemma 3.3, and since uj is bounded away from zero 
to the past by Lemma 7.4. Assume uj does not converge to ui. Then there is a 
sequence Sk —00 such that uj{sk) — + uj2 where we can assume uj2 > uji. We can 
also assume < Sk- Then 

uj{sk) = exp( / {q- 2)ds)u){Tk). 

Since 

q-2< 3{NiN2 + N2N3 + N3N1) < 3e"^ 



for T < T by Lemma 7.5 and the constraint (2.3), we have, assuming Sk < T, 

3 



Thus 



/ {q- 2)ds < 3 / e°"'dT < 

J Ti. J Ti. 



uj{sk) < exp(-e"^''-)w(rfc) uJi, 
a 



so that UJ2 < uji contradicting our assumption. □ 

Corollary 7.2. Consider a solution to ^2^)-^^) with 7 = 2 and uj > 0. Then 

lim N,{t) = 

/ori = 1,2,3. 

Proof. Assume iVi does not converge to zero. Then there is a type II a-limit point 
with iVi and uj non-zero by Corollary 7.1 and Lemma |7.6| . If we apply the flow, we 
get a-limit points with different uj in contradiction to Lemma [7.(^. □ 



Lemma 7.7. Consider a solution to (2.1)-(2.S) with 7 = 2 and a; > 0. // it has 



an a-limit point of type I inside the triangle, the solution converges to that point. 

Proof. Let x be the limit point. Let _B be a ball of radius e in wl]-|_l]_-space, with 
center given by the w, S_|_, E_-coordinates of x. Let Tk — > —00 be a sequence that 
yields x. Assume the solution leaves B to the past of every Tk. Then there is a 
sequence Sk — *■ —00, such that the uj, S+, E_ -coordinates of the solution evaluated 
in Sk converges to a point on the boundary of B, Sk < Tk, and the E-|_,E_- 
coordinates of the solution are contained in B during [sk,Tk], k large enough. 
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Since all expressions in the Ni decay exponentially as e"'^, for some a > 0, as long 
as the Lo, I]_ -coordinates are in B (e small enough), we have 



IS'. 



for r G [sfcjTfc] where at — > 0. We get 



|E+(rfe)-E+(sfc)|<^ 
a 



0, 



and similarly for S_ and lo. The assumption that we always leave B consequently 
yields a contradiction. We must thus converge to the given a-limit point. □ 



Proposition 7.1. Consider a solution to (2.1 )-{2.c ) with ^ — 2 and uo > Q. If Ni 

is non- zero for the solution, it converges to a type I point in the complement of Mi 
with a; > 0. 



Proof. If there is an a-limit point on Mi, we can use Lemma 7.2 to obtain a 
contradiction to Lemma 7.6. If there is an a-limit point in Mk and is zero 
for the solution, the solution converges to that point by an argument similar to 
the one given in the previous lemma. What remains is the possibility that all the 
a-limit points are on the Ck. Since uj converges, the possible points projected to the 
I]+I]„-plane are the intersection between a triangle and a circle. Since the a-limit 
set is connected, we conclude that the solution must converge to a point on one of 
the Ck- □ 



Proposition 7.2. Consider a solution to ( \2.i )- (i2.3i ) with 7 = 2 and uj > 0. If Ni 

is non-zero for the solution, the solution cannot converge to a point in Ci . 



Proof. Assume i — \. Then Ci is the subset of (7^) consisting of points with 
S+ = 1/2 and > 0. Since N2, N3, N2N3, N2N1 and N3N1 converge to zero 
faster than N^, will in the end be positive, cf. (7^), so that there is a T such 
that S_|_(r) > 1/2 for T < r. Since A^i will dominate in the end, we can also assume 
^(t) < 2 for T <T. By (2.1) we conclude that jA^il increases backward as r < T 
contradicting Corollary 7.2. □ 

Adding up the last two propositions, we conclude that the S-|_S_-variables of 
Bianchi VIII and IX solutions converge to a point interio r to the triangle of Figure 
H, and Q to the value then determined by the constraint (|2.3| ). In the Bianchi VIIq 
case, a side of the triangle disappears, increasing the set of points to which S+, S_ 
may converge. We sum up the conclusions in Section |l^. 



8. Type I solutions 

Consider type I solutions {Ni = 0) . The point F and the points on the Kasner circle 
are fixed points. Consider a solution with < ^{tq) < 1. Using the constraint, we 
may express the time derivative of Q in terms of Q. Solving the resulting equation 
yields 

lim n{T) = 0, lim r2(r) = 1. 



By (|2T|) (S+,S„) moves radially. 
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Proposition 8.1. For a type I solution, with 2/3 < 7 < 2, which is not F, we have 
hm (E+,I]_,r!)(T) = (a+/|a|,a_/|a|,0), 



r — ^ — 00 



where (ct+,(7_) is the initial value of and \a\ is the Euclidean norm of 

the initial value. 



9. Type II solutions 

Proposition 9.1. Consider a type II solution with Ni > and 2/3 < 7 < 2. // 

the initial value for E_ is non-zero, the a-limit set is a point in IC2 U /C3. // the 
initial value for E_ is zero, either the solution is the special point P^[II), it is 
contained in Tu, or 



(9.1) lim (f],E+,iVi)(r) = (0,-1,0). 

r — > — oo 

Proof. Let the initial data be given by (it+, f7_, fJo)- Tiie vacuum case was liandled 



in Proposition 5.1, so we will assume Q,q > 0. 



Consider first the case (T_ 7^ 0. Compute 

Thus, I]_ decreases if it is negative, and increases if it is positive, as we go back- 



ward in time, by (2.1). Thus, both A^i and must converge to as t ^ —00, 
since the variables are constrained to belong to a compact set, and because of the 
monotonicity principle. Since E_ is monotonous and the a-limit set is connected. 



see Lemma 3.1, (E+,E_) must converge to a point, say (s+,s_) on the Kasner 



circle. We must have s_ 7^ 0, and 

2sl_ + 2s^ - 4s+ > 0, 

since iVi converges to 0. There are two special points in this set, but we may 
not converge to them, since that would imply Ni = for the entire solution by 



Proposition 3.1. The first part of the proposition follows. 



Consider the case cr_ = 0. There is a fixed point P^{II). Eliminating from 



(2.1)-(2.3), we are left with the two variables iVi and E+. The linearization has 
negative eigenvalues at P^{II), so that no solution which does not equal P^{II) 
can have it as an a-limit point, cf. |jl^ pp. 228-234. There is also a set of solutions 
converging to the fixed point F. Consider now the complement of the above. The 
function 



(l-«S+)2' 

where v = (37 — 2)/8 and m — 3w(2 — 7)/8(l — u^), found by Uggla satisfies 

3(2-7) 1 2 
^^-l-.E+l-.2(^+-^) 

Apply the monotonicity principle. Let G — and U be defined as the subset of 
riE+TVi-space consisting of points different from P^ (IT), which have > 0, A^i > 
and < 1. Let M be defined by the constraint. If S+ = v then Z7 = 0, but if 
we are not at P^{II), E4. = v implies S'^. ^ 0. Thus, G o x is strictly monotone 
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as long as x is contained in [/ n M . Since the solution cannot have P^{II) as an 
a-limit point, we must thus have A^i = or — in the a-limit set. Observe that 

(9.2) SV = ^7V2(2_s+)-^(2-7)J7I]+. 

Thus, if the solution attains a point < 0, then holds. We will now prove 
that this is the only possibility. 

a. Assume we have an a-limit point with Ni > and = 0. Then we may apply 
the flow to that limit point to get = — 1 as a limit point, but then the solution 
must attain < 0. 

b. If > but iVi — 0, then we may assume S+ ^ since we are not on jTu, cf. 



Lemma 4.2. Apply the flow to arrive at S_|_ = — 1 or !]_(_ = 1. The former alternative 
has been dealt with, and the latter case allows us to construct an a-limit point with 
A^i > and = 0, since iVi increases exponentially, and decreases exponentially, 
in a neighbourhood of the point on the Kasner circle with S-|_ = 1, cf. Proposition 



3.1 



c. The situation il = A^i = can be handled as above. □ 

We make one more observation that will be relevant in analyzing the regularity of 
^11. 

Lemma 9.1. The closure of Tu does not intersect A. 

Proof. Assume there is a sequence Xk G J'li such that the distance from to A 
goes to zero. We can assume that all the x^. have A^i > by choosing a suitable 
subsequence and then applying the symmetries. We can also assume that Xk x G 



A. Since S]_ = for all the Xk by Proposition 3.1, the same holds for x. Observe 
that no element of can have < 0, because of (^^). If Ni corresponding to 
X is zero, we then conclude that x is defined by S+ — 1 and all the other variables 
zero. Applying the flow to the past to the points Xk will then yield a sequence 
Tjk G ^11 such that i/k converges to a type II vacuum point with A'^i > and 



E_ = 0, cf. the proof of Proposition 3.1. Thus, we can assume that the limit point 



X € A has iVi > 0. Applying the flow to x yields the point S+ = — 1 on the Kasner 



circle by Proposition 5.1. By the continuity of the flow, we can apply the flow to 



Xk to obtain elements in J-u with S_|_ < which is impossible. □ 

10. Type VIIq solutions 

When speaking of Bianchi VIIq solutions, we will always assume A'^i = and 
N2,Ns> 0. Consider flrst the case iV2 = N3 and !]_ = 

Proposition 10.1. Consider a type VIIq solution with Ni = and 2/3 < 7 < 2. 

If N2 — Nt, and S_ = 0, one of the following possibilities occurs 

1. The solution converges to ~ 1 on the Kasner circle. 

2. The solution converges to F. 

3. limT-^_oo = —1, liniT-^-oo -^2 — n2 > 0, limT-^_oo = 0. 



Proof. Since 



S'+ = -^(2-7)f^SH 
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if N2 — N3, the conclusions of the lemma follow, except for the statement that N2 
converges to a non-zero value if S+ converges to —1. However, fl will decay to zero 
exponentially close to the Kasner circle, and by the constraint, 1 + E+ will behave 
as n close to E-)_ = — 1. Thus, q + 2S_|. will be integrable. □ 

Before we state a proposition concerning the behaviour of generic Bianchi VIIo 
solutions, let us give an intuitive picture. Figure ^ shows a simulation with 7=1, 
where the plus sign represents the starting point, and the star the end point, going 
backward. Q will decay to zero quite rapidly, and the same holds for the product 
N2N3. In that sense, the solution will asymptotically behave like a sequence of 
type II vacuum orbits. If both N2 and N3 are small, and we are close to the section 
IC2 on the Kasner circle, then A^2 will increase exponentially, and N3 will decay 
exponentially, yielding in the end roughly a type II orbit with N2 > 0. If this orbit 
ends in at a point in /C3, then the game begins anew, and we get roughly a type II 
orbit with > 0. Observe however that if we get close to /Ci, there is nothing to 
make us bounce away, since A'^i is zero. The simulation illustrates this behaviour. 
Consider the figure of the solution projected to the E+E_-plane. The three points 
that appear to be on the Kasner circle are close to JC2, IC3 and JCi respectively. 
Observe how this correlates with the graphs of N2, N3 and q. 




-2 -1 1 2 5 10 15 




5 10 15 5 10 15 



-T -X 

Figure 4. Illustration of a Bianchi VIIq solution. 



Proposition 10.2. Generic Bianchi VIIq solutions with Ni — and 2/3 < 7 < 2 

converge to a point in K-i. 
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We divide the proof into lemmas, 
contained in a compact set. 



First we prove that the past dynamics are 



Lemma 10.1. For a generic Bianchi VIIq solution with Ni = and 2/3 < j < 2, 

(A^2, -^3)(— oo, 0] is contained in a compact set. 



Proof. For a generic solution, 



Z-i 



is never zero. Compute 
(10.1) 



16 



N2N3 
S^(l + S+ 



3 ^tv2 



:Z-i. 



gtSi + (iV2 - N3)- 

The proof that the past dynamics are contained in a compact set is as in Rendall 
ITl. Let r < 0. Then 



so that 



Z_i(r) >Z_i(0), 
(^2^3) (t) < 



3^-1 (0) 

Combining this fact with the constraint, we see that all the variables are contained 
in a compact set during (— cx), 0]. □ 

We now prove that N2N3 0. The reason being the desire to reduce the problem 
by proving that all the limit points are of type I or II, and then use our knowledge 
about what happens when we apply the flow to such points. 

Lemma 10.2. Generic Bianchi VIIq solutions with Ni = and 2/3 < 7 < 2 
satisfy 

lim (7V2iV3)(T) =0. 



Proof. Assume the contrary. Then we can use Lemma 10.1 to construct an a-limit 



point (cj, (T+, (T_, 0, rt2, 'T's) where 71.2^3 > 0. We apply the monotonicity principle 
in order to arrive at a contradiction. With notation as in Lemma 5.1, let U be 
defined by N2 > 0, N3 > and + - Ns)^ > 0. Let G be defined by 
and M by the constraint ( 2_^ ) . We have to show that G evaluated on a solu tion i s 
strictly monotone a s lo ng as the solution is contained in U O M. Consider ( 10.1 ). 
By the constraint (^!3|) , + {N2 - iVa)^ > imphes T.+ > -1. Furthermore, 
Z_i > on C/. If Z'_i = in [/ n A/, we thus have S- = 0, but then S'_ ^ since 
^- + {N2 - N^)'^ > and A^2 + ^3 > 0. The a-limit point we have constructed 
cannot belong to U. On the other hand, 712, > and since Z_i increases as we 
go backward, + (^2 — ".3)^ cannot be zero. We have a contradiction. □ 



Proof of Proposition 10. L Compute 
(10.2) = -(2- 217- 2S^ - 



2E^)(1 + E+) - -(2-7)^7E^ 



by (2.4). Assume we are not on 'Pviio oi J-yug. Let us first prove that there is an 
a-limit point on the Kasner circle. Assume F is an a-limit point. Then we may 
construct a type I limit point which is not F, and thus a limit point on the Kasner 



circle, cf. Lemma 4.2 and Proposition 8.1. By Lemma 10.2, we may then assume 
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that there is a hmit point of type I or II, which is not P^{II) or P^{II), and does 
not he in J-'i or J-u, cf. Lemma 4.1. Thus, we get a hmit point on the Kasner circle 



by Proposition 8.1 and Proposition 9.1 



Next, we prove that there has to be an a-hmit point which hes in the closure of /Ci. 
If the a-limit point we have constructed is in IC2 or /C3, we can apply the Kasner 
map according to the remark following Proposition |6.l| . After a finite number of 
Kasner iterates we will end up in the desired set. If the a- limit point we obtained 
has S+ = —1, we may construct a limit point with 1 + E+ = e > by Proposition 
3.1. We can also assume that J7 = for this point, since decays exponentially 
going backward when is close to —1. By Lemma 10.2, this limit point will be 
a type I or II vacuum point, and by applying the flow we get a non special limit 
point on the Kasner circle. As above, we then get an a-limit point in the desired 
set. Let the E+S_-variables of one a-limit point in the closure of JCi be (cr_|_,cr_). 

By (10.2), we conclude that once S+ has become greater than 0, it becomes mono- 
tone so that it has to converge. Moreover, we see by the same equation that fl then 
has to converge to zero, and + has to converge to 1. Since the a-limit set is 
connected, by Lemma [3.l| and Lemma 10.1, we conclude that S_) has to con- 
verge to ((7+, (J-). By Proposition 3^, ((t+,(t_) cannot equal (1/2, ±\/3/2), since 
otherwise N2 or A^3 would be zero for the entire solution. Consequently, (T+ > 1/2, 
and we conclude that N2 and A^3 have to converge to zero. The proposition follows. 
□ 



11. TAUB type IX SOLUTIONS 

Consider the Taub type solutions: E_ = and N2 — N3. We prove that except for 
the cases when the solution belongs to JFjx or Vix, (S+, E^) converges to (—1,0). 

Lemma 11.1. Consider a type IX solution with E_ — 0, N2 — N3 and 2/3 < 7 < 
2. Then I]_|_(to) < and 57 (tq) < 1 imply 

lim (r!,E+,E_,7Vi,iV2,7V3)(T) = (0, -1, 0, 0, na, na), 
where < ?i2 < 00. 

Proof. We prove that the flow will take us to the boundary of the parabola 51 -I- E^ = 
1 with E_|_ < 0, and that we will then slide down the side on the outside to reach 
S+ = — 1, see Figure |. The plus sign in the figure represents the starting point, 
and the star the end point. 

1. Let us first assume E+(ro) < 0, 51(ro) < 1 and 51(to) -I- E^(to) > 1. Consider 

c = {T<To:te[T, To] ^ E+(t) < 0, n{t) < n{To), n{t) + T^Kt) > i}. 

We prove that C is not bounded from below. Assume the contrary. Let t be the 
infimum of C, which exists since C is non-empty and bounded from below. Since 
teC, E+(i) < 0. Let t' <the such that E+ < in [t',t]. Observe that 

(11.1) n' = [(37 - 2)(f7 + E^ - 1) + 3(2 - 7)E^]f7. 
By the constraint, 

(11.2) n + ^l-i = ^NU^^-i). 
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Figure 5. Part of a Taub type IX solution projected to the E+fi-plane. 



Since < in [t' , t], N2/N1 increases as we go backward in that interval, because 
of 



Consequently + > 1 in [t' , t], by ( |ll.2 ), so that decreases in the interval by 
(11.1). Thus t' E C, contradicting the fact that t is the infimum of C. 



No 



Let r < To. Then S+(r) < -^1 - n{To). By ( |ll.lD , wc then conclude n ^ 0. By 
(PH), we also conclude that N1N2 -> and Ni 0. By ( |ll.2[ ), we have Y,+ -1. 
Using the constraint (11.2) and ( |2.2| ), we conclude that q + 2S+ is integrable, so 
that N2 = N3 will converge to a finite non-zero value. 

2. Assume now S+(to) < 0, ri(To) < 1 and ^{tq) + S^(to) < 1. Observe that 



(11.3) 



= (1 - O - S^)(4 - 2E+) - -(2 - 7)r!S+ + 9N1N2 



As long as + < 1, E+ decreases as we go backward in time by ( 11.3 ). Then 



N2/N1 will increase exponentially until il + E^ = 1, by the constraint, and < 0. 

□ 

Lemma 11.2. Consider a type IX solution with E_ —Q,N2 — N3 and 2/3 < 7 < 
2. It is contained in a compact set for r < and N1N2 ^ 0. 



Proof. Note that iVi must be bounded for r < 0, as follows from Lemma |3.3| , the 
fact that N2 — N3, and the fact that N1N2N3 decreases backward in time. To 
prove the first statement, assume the contrary. Then there is a sequence Tk —00 
such that N2{Tk) 00. We can assume N2{Tk) < 0, and thus 



(11.4) 



i(37-2)f7 + 2S^ + 2E+ <0 
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in Tfc. Since N1N2 is decreasing as we go backward, A^i and N1N2 evaluated at Tk 
must go to zero. Thus ft + — 1 will become arbitrarily small in by (11.2). If 
ft{Tk) > 1 for all k, we get 

S+(r,)<-i(37-2) 



by (11.4), so that 



which is a contradiction. In other words, there is a /c such that E-i-(rfc) < 0, by 



(11.4), and fl{Tk) < 1. We can then use Lemma 11.1 to arrive at a contradiction 
to the assumption that the solution is not contained in a compact set. 

To prove the second part of the lemma, observe that iViiV| converges to zero, as 



follows from the existence of an a- limit point and Lemma 5.2. Thus 
N1N2 = iVy^[7Vi7V|]i/2 < C[iViiV|]i/2 _^ 0. 

□ 

Proposition 11.1. For a type IX solution with S_ = 0, A^2 = and 2/3 < 7 < 

2, either the solution is contained in J-'ix or Vix, or 



lim {n,j:+,j:^,Ni,N2,N3){T) = (0,-1,0,0,712,712) 

r — > — 00 

where < rt2 < oo. 

Remark. Compare with Proposition |3.l| . Observe also that when S+ for the solution 
converges to —1, we approach S+ ~ — 1 , il — from outside the parabola il + E^ = 
1, as follows from the proof of Lemma [l 1 . l| . 

Proof. Consider a solution which is not contained in J-ix or Vix- By Lemma |ll.2 



there is an a- limit point with N1N2 = 0. We can assume it is not P^{II). We 
have the following possibilities. 

1. It is contained in JFj U J^n U JViIq ■ Then F is an a-limit point. Since the solution 



is not contained in Tix, we get a type I limit point which is not F, b y L emma 4.2 
and thus either = — 1 or = 1 as limit points, by P ropo sition ^.l| . The first 
alternative implies convergence to — —1, by Lemma |ll.l| . If we have a ty pe 



I a-limit point with = 1, we can apply the Kasner map by Proposition 3.1 
order to obtain a type I limit point with = — 1. 

2. The limit point is of type I. This possibility can be dealt with as above. 

3. It is of type II. We can assume that it is not P^{IT), by Lemma 4.1, and that 



it is not contained in jFjj. Thus we get = — 1 on the Kasner circle as an a-limit 



point, by Proposition 9.1, and thus as above convergence to = — 1. 



4. The limit point is of type VIIq. We can assume 7^ 0. If < 0, we can 
apply Lemma [l 1 . l| again, and if S+ > 0, we get = 1 on the Kasner circle as an 



a-limit point, by Proposition 10.1, a case which can be dealt with as above. □ 
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12. Oscillatory behaviour 



It will be necessary to consider Bianchi IX solutions to (2.1)-( |2.3| ) under circum- 
stances such that the behaviour is oscillatory. This section provides the technical 
tools needed. 



Let g he a, function, 

(12.1) A = 

and X = (i, y)* satisfy 



9 
'9 



X = Ax + e, 
where e is some vector valued function. 

Lemma 12.1. Let (j)Q be such that (sin(0o)i cos(0o)) '^'^'^ {S;(TQ),y{TQ)) are parallel. 
Define 

(12.2) C(t) = / 9is)ds + ^o 

J To 

and 

(12.3) x(r) = 



x{t) \ ^ ( sinC(r) 
y(r) ) I cos.^(r) 



Then 

(12.4) ||x(r)-x(r)|l<|l-(.^2(ro)+y2(ro))i/2| + | / \\e{s)\\ds\. 



Proof. Let 



X y 

We have [A, <I>] = 0, <!>' = -A^ and x' = Ax. We get 

($(x - x))' = -A$(i - x) + $(A(A - x) + e) = $e. 

Thus 

(x - x)(t) - $-i(T)$(ro)(i - x)(ro) + ^-^(t) / $(s)e(s)ds 



But <i> takes values in 5(9(2) and the lemma follows. □ 

In order to prove the existence of an a-limit point for Bianchi IX solutions, and 
that, generically, there is a limit point on the Kasner circle, we need the following 
lemma. 

Lemma 12.2. Consider a Bianchi IX solution with 2/3 < 7 < 2. Assume there is 
a sequence —00 such that qijk) 0, and N2(Tk), N^^Tk) 00, then for each 
T , there is a t < T such that S+(t) > 0. 



Proof Observe that by (|2J), q^O and N2 + N3> Ni implies < -2. However, 
the only term appearing in the constraint which does not go to zero in Tk is {N2 — 
Ns)"^, since the product N1N2N3 decreases as we go backward. Thus |I]'_(Tfe)| 00, 
and the behaviour is oscillatory. It is clear that could become positive during 
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the oscillations, but only when |S]_| is big, so that we on the whole should move in 
the positive direction. 

Assume there is a T such that '^+{t) < for all t <T. 
We begin by examining the behaviour of different expressions in the sets 

=U^^,[r„-l,T„] 

and 

2? = U,r=i[r„-l,T„]. 

Observe that by the fact that {il, S+, £ „) a re constrained to belong to a compact 
set during (—00, 0], according to Lemma 3.3, N2 and N3 go to infinity uniformly in 
V (by which we will mean the following): 

VAf 3K ■.k>K ^ Niir) > M Vr e Vk, i = 2, 3. 



Thus iVi and Ni{N2 + N3) go to zero uniformly in V. By (2J), Q also converges to 
zero uniformly i n V . Due to the constraint, we get a bound on + |(A^2 — -^3)^ 
in D. Consider ( ^.4| ). The last two terms go to zero uniformly. If the first term is 
not negative, 1 — f2 — — < 0. By the constraint, it will then be bounded by 
an expression that converges to zero uniformly in V. Thus, for every 6 > there 
is a such that k > K implies < S in Vk- Combining this with the fact that 
Q{Tk) ^ 0, and the assumption that S+(r) < for t < T, we conclude that S]-|_ 
converges uniformly to zero in V. 



Next, we use Lemma 12.1 in order to approximate the oscillatory behaviour. Define 
the functions 



y 



V3 N2 



1/2 
-N3 



We can apply Lemma 12.1 with 
g = -3{N2 + 



2 (1 



2(1 



+ = gi +92 



and ex, f-y given by (15.5) and (15.6), cf. Lemma 15.1. By the above, we conclude 
that X and y are uniformly bounded on P^, if k is great enough, and that ||e| 



converges to zero uniformly on D. Let be the expression given by Lemma 12.1 
with To replaced by tj, and <po by a suitable (pk- Let S > 0. By the above and 
q{Tk) 0, we get 

(12.5) ||(i-xfe)(r)|| <5, 

if T G [Tfc — 1, Tfe], and k is great enough. In [tj. — 1, r^], we thus have 

(12.6) SV = -2 + 24(l-S^) + Pfc, 

where the error pk can be assumed to be arbitrarily small by choosing k great 
enough, cf. (2.4). 

Let 

ffc(T) = / g{s)ds + 4>k 



THE BIANCHI IX ATTRACTOR 



31 



be as in (12.2). Since N2 + N3 goes to infinity uniformly, [r^ — 1, r^] can be assumed 
to contain an arbitrary number of periods of ^fe, if A: is great enough. Thus, we 
can assume the existence of ri^fe,r2,fe € [Tfc — l,Tfe], such that T2,fe — Ti^k > 1/2 
and ^k{Ti.k) - £,k{T2,k) is an integer multiple of tt. Let [ri,T2] C [Ti^k,T2,k] satisfy 
^k{Ti) — ^k{T2) = TT. We can assume T2 — n to be arbitrarily small by choosing k 



great enough. Considering (2.1), and using the fact that q is bounded, we conclude 
that N2 + N3 cannot change by more than a factor arbitrarily close to one during 
[''"ii ''2]- Since the expression involving iV2 + A'3 dominates g, we conclude that 

3 

-^.9(Tmax) < -g(Tmin), 

where T^ax and Tmin correspond to the maximum and the minimum of — g in [ti, T2]. 
Estimate 

^■"0^2., ^2^.„ 24(1 -s^)^^ 2x1(1 



/ 2x1(1 -T.l)ds^ ±ldr] = - 

< -/ 2sin^(i])df] ^ -. 

g(Tmin) JCfc(ri)-7r 9{'''min) 

We get 

/•?fc(-r2) 2 TT 3 /""^^ 

T2 - n = / -dr; > > - / 2x1(1 - ^l)ds. 

Jikiri) 9 5(Tmax) 4 

Consequently, ( 12. 6| ) yields 



S+(t2) - S+(ti) = -2(t2 - Ti) + / 24(1 - E^)ds + /" pkdT < 



T2 



< "^('^2 - n) + / Pfcrfr. 



Since ^fc(Ti.fc) — ^fe(T2^fc) corresponds to an integer multiple of tt, we conclude that 

2 r'''2.k 2 /■■''2,fc 

S+(T2,fe) - i;+(Ti_fe) < --(r2,A; - Ti.fc) + / Rkd-T < - - + I pkdr. 

However, the expressions on the far left can be assumed to be arbitrarily small, and 
the integral of pk can be assumed to be arbitrarily small. We have a contradiction. 
□ 



13. BlANCHI IX SOLUTIONS 

We first prove that there is an a-limit point. If we assume that there is no a-limit 
point, we get the conclusion that the EucHdean norm || A^|| of the vector (Ni, N2, N3) 
has to converge to infinity, since (fl, E+, E_) is constrained to belong to a compact 
set to the past by Lemma [3.3| . In fact. Lemma ^.3| yields more; it implies that two 
Ni have to be large at any given time. Since the product N1N2N3 decays as we go 
backward, the third Ni has to be small. Sooner or later, the two Ni which are large 
and the one which is small have to be fixed, since a 'changing of roles' would require 



two Ni to be small, and thereby also the third by Lemma 3.3, contradicting the 
fact that ||A^|| — > 00. Therefore, one can assume that two Ni converge to infinity, 
and that the third converges to zero. More precisely we have. 
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Lemma 13.1. Consider a Bianchi IX solution. //||A'^|| — > oo, we can, by applying 
the symmetries to the equations, assume that N2, -/V3 — s- 00 and Ni, Ni{N2+N3) — s- 
0. 

Proof. As in the vacuum case, see iQ. □ 

Lemma 13.2. A Bianchi IX solution with 2/3 < 7 < 2 has an a-limit point. 

Proof. If the solution is of Taub type, we already know that it is true so assume 
not. We assume N2, — s- 00, since if this does not occur, there is an a-limit 



point by Lemma 3.3 and Lemma 13.1. By ( p.4[ ) we have < if S+ = using 
the constraint (assuming A^2 + -^3 > 3A^i). Thus, there is a T such that if E+ 
attains zero in r < T, it will be non-negative to the past, and thus N2N3 will be 
bounded to the past since S+ has to be negative for the product to grow . If there 



is a sequence —00 such that g(rfc) 0, we can apply Lemma 

at a contradiction. Thus there is an S such that 



12.2 to arrive 



(13.1) 

for all T <S. 
Consider 

(13.2) 



Z-i 



q{T) >€>{) 



+ {N2 



N2N3 

The reason we consider this function is that the derivative is in a sense almost 
negative, so that it almost increases as we go backward. On the other hand, it 
converges to zero as r — > —00 by our assumptions. The lemma follows from the 
resulting contradiction. We have 



(13.3) 

Letting 



f E^(l 



4:V3^-{N2 ~ N3)Ni 



N2N3 



N2N. 



/=-I]2_ + (7V2-iV3)', 



we have, using the constraint, 

h < 4T.^_Ni{N2 + N3) + 2V3Nif < NiN2N3f 
for, say, t <T' < S. Thus 

(13.4) Z'_^ < NiN2N3Z^i 



for all T < T'. Since q > e > for aU r < T' < S" by ( |l3.l| ), we get 

(iViiV2iV3)(r) < {NiN2N3){T')exp[3e{T-T')] 



for T < T' . Inserting this inequality in ( 13.4 ) we can integrate to obtain 



Z-i{r) > Z_i(r')cxp(--(iVi7V2iV3)(T')) > 

3e 



for T < T'. But Z_i(r) 
contradiction. □ 



as T 



-00 by our assumption, and we have a 
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Corollary 13.1. Consider a Bianchi IX solution with 2/3 < 7 < 2. For all e > 0, 

there is a T such that 



for all T <T. Furthermore 



17 + + S?. < 1 + e 



lim iNiN2N3){T) = 0. 

T — * — 00 



Proof. As in the vacuum case, see |16 . The second part follows from Lemma 
and Lemma 13.2. □ 



Proposition 13.1. A generic Bianchi IX solution with 2/3 < 7 < 2 has an a-limit 
point on the Kasner circle. 



Proof. Observe that by Lemma 13.2 and Corollary 13.1, there is an a-limit point 
of type I, II or VIIq. 

1. First we prove that we can assume the a-limit point to be a type VIIq point with 
iVi = 0, < iVa = A^3, n^O, E_ ^ and S+ = -1. 

a. If there is an a-limit point in J-j, J-u or J-yu^, F is a limit point, but then there 



is an a-limit point on the Kasner circle, by Lemma L2 and Proposition 8.1 . 

b. Assume there is an a-limit point in Vyug, or that one of P^{II) is an a-limit 



point. Then there is a limit point of type II which is not i-y (//), by Lemma 4.1 



and we can assume it does not belong to Tn. We thus get an a-limit point on the 



Kasner circle by Proposition 9.1 



c. Consider the complement of the above. We have an a-limit point of type I, II 
or VIIq which is generic or possibly of Taub type. If the limit point is of type I or 
II, we get an a-limit point on the Kasner circle by Proposition |8.l| and Proposition 
|9.l| . If the limit point is a non-Taub type VIIq point, we get an a-limit point on 
the Kasner circle by Proposition 10.2 . Assume it is of Taub type with I]_ = 0, 
A^2 = N^. By Proposition 10.1, we can assume that we have an a-limit point of the 
type mentioned. 

2. We construct an a-limit point on the Kasner circle given an a-limit point as in 
1. Since the solution is not of Taub type, we must leave a neighbourhood of the 
point £_) = (—1, 0). If A^2 and TVs evaluated at the times we leave do not go 
to infinity, we are done. The reason is that we can choose the neighbourhood to be 
so small that and iVi decrease exponentially in it, see (2.1). If N2{tk) or N^ltk) 
is bounded, we get a vacuum Bianchi VIIq a-limit point which is not of Taub-type 
by choosing a suitable subsequence (if we get a type I or II point we are done, see 
the above arguments). By Proposition 10.2, we then get an a-limit point on the 
Kasner circle. Thus, we can assume the existence of a sequence tk — > —00 such that 
N2{tk) and N3{tk) go to infinity. 

There are two problems we have to confront. First of all iV2 and have to decay 
from their values in in order for us to get an a-limit point. Secondly, and more 
importantly, we need to see to it that we do not get an a-limit point of the same 
type we started with. Let us divide the situation into two cases. 



34 



HANS RINGSTROM 



a. Assume that for each tk there is an Sk < tk such that I]+(sfe) = 0. Observe that 
when = 0, we have 

K < Ini{9Ni - SN2 - 3N3) 



by the constraint ( |2.3| ), and {2A). Thus, we can assume that we have 3iVi > N2+N3 
in Sk, since there is an a-hmit point with S+ = — 1. Thus there must be an < t/c 
such that, at r^, either Ni = N2 < N3, Ni ^ N3 < N2 or iVi < 7V2, A^i < N3 and 
3A^i > N2 + N3. One of these possibihties must occur an infinite number of times. 
The first two possibihties yield a type I or II hmit point, and the last a type I limit 
point because, of the fact that N1N2N3 and Lemma 3.3. As above, we get an 
a-limit point on the Kasner circle. 

b. Assume there is a T such that S+(t) < for all t < T. Then iVi 0, since 
Ni{tk) 0, and S+ < implies that A^i is monotone. Assume there is a sequence 
Tk — !■ —00 such that N2 or N3 evaluated at it goes to zero. Then we get an a limit 
point of type I or II, a situation we may deal with as above. Thus we may assume 
Ni > e > 0, i = 2,3 to the past of T. Similarly to the proof of the existence of an 
a-limit point, we have 

Z'_i < c,NiN2N3Z^i. 

If there is an S and a ^ > such tha t ^(t ) > ^ > for all t < S", we get a 
contradiction as in the proof of Lemma 13.2 , since (A^2A^3)(ife) ~^ 00 ■ Thus there 
exists a sequence rj. —00 such that q(Tk) 0. If N2(Tk) or N3{Tk) contains a 
bounded subsequence, we may refer to possibilities already handled. By Lemma 



12.2, we get > 0, a contradiction. □ 



14. Control over the density parameter 



The idea behind the main argument is to use the existence of an a-limit point on 
the Kasner circle to obtain a contradiction to the assumption that the solution does 
not converge to the closure of the set of vacuum type II points. The function 

d=n + N1N2 + N2N3 + N3N1 

is a measure of the distance from the attractor. We can consider d to be a function 
of T, if we evaluate it at a generic Bianchi IX solution. If Tk —00 yields the 
a-limit point on the Kasner circle, then d{Tk) —^0. If d does not converge to zero, 
then it must grow from an arbitrarily small value up to some fixed number, say 
(5 > 0, as we go backward. In the contradiction argument, it is convenient to know 
that the growth occurs only in the sum of products of the Ni, and that during 
the growth one can assume fl to be arbitrarily small. The following proposition 
achieves this goal, assuming S is small enough, which is not a restriction. The proof 
is to be found at the end of this section. 

Proposition 14.1. Consider a Bianchi IX solution with 2/3 < 7 < 2. There exists 
an e > such that if 

(14.1) N1N2 + N2N3 + N1N3 < e 

in [ti, T2], then 

n < c^n{T2) 

in [ti,T2] if fl(T2) < e. Here c-y > only depends on 7. 
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The idea of the proof is the following. If the sum of product of the Ni and fl are 
small, the solution should behave in the following way. If all the Ni are small, 
then we are close to the Kasner circle and decays exponentially. One of the Ni 
may become large alone, and then increases, but it can only be large for a short 
period of time. After that it must decay until some other Ni becomes large. But 
this process of the Ni changing roles takes a long time, and most of it occurs close 
to the Kasner circle, where fl decays exponentially. Thus, may increase by a 
certain factor, but after that it must decay by a larger factor until it can increase 
again, hence the result. Figure illustrates the behaviour. 




-2 -1 1 2 0.5 1 1.5 2 2.5 




0.5 1 1.5 2 2.5 0.5 1 1.5 2 2.5 



-T -T 

Figure 6. Part of a type IX solution. 

We divide the proof into lemmas, and begin by making the statement that Q decays 
exponentially close to the Kasner circle more precise. 

Lemma 14.1. Consider a Bianchi IX solution with 2/3 < 7 < 2. // 

^l + ^l > ^(37 + 2) 

in an interval [si,S2], then 
for s e [si, 52]; where 

«7 = ^(2-7)- 
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Proof. Observe that 
(14.2) 



f7'>4[S^ + S2_--(37-2)]0, 



so that under the conditions of the lemma 



The conclusion follows. □ 



Next , we prove that if the Ni all stay sufficiently small under a condition as in 
( L4.1 ) and Q starts out small, then Q will remain small. 

Lemma 14.2. Consider a Bianchi IX solution with 2/3 < 7 < 2. There is an 
e > such that if 



|Ar2< i(6_3^) 



(14.3) 

(14.4) N1N2 + N2N3 + N1N3 < e 

in an interval [si,S2\, and f2(s2) < then Q(s) < fl{s2) for all s G [si,S2]. 
Proof. Let 

£ = {r e [si, S2] : t € [r, S2] ^ n{t) < n{s2)}. 

Let T £ £, T > si. There must be two Ni, say N2 and i Va, su ch that A''2 < e-^^^ and 
N3 < e^/^ in r, by (|l4.4|) . By the constraint ( [2!^ ) and ( |l4.4| ), we have in r, 

+ > 1 - -A^2 _ f7 _ > 1(3^ + 2) - 4e, 
4 8 



so that assuming e small enough depending only on 7, we have ^'{t) > 0, cf. (14.2). 
Thus there exists an s < r such that s e £. In other words, £ is an open, closed, 
and non-empty subset of [si, S2], so that £ = [si, S2]. □ 

The next lemma describes the phase during which Q may increase. 

Lemma 14.3. Consider a Bianchi IX solution with 2/3 < 7 < 2. There is an 
e > such that if 



(14.5) 
(14.6) 



3^r2> i(6_37) 
iViiV2 + N2N3 + N1N3 < e 



in [si,S2], and il(s2) < e, then S2 — Si < Ci,j and Q{s) < C2,-yil(s2) for all s € 
[31,82], where Ci.^ and 02.^ are positive constants depending on 7. 

Proof. Assume e is small enough that 



so that A^i > e^/^ in [si, 52]- Assuming e < 1 we get Ni < e^l"^ in \s\, S2], i = 2, 3. 
Use the constraint (2^) to write 



(14.7) 



where < 3e by ( |14.6| ). Thus, 

1-17- 



> ^^1/2 _ 
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SO that we may assume 

(14.8) fi + S^ + S?.<l 

in [si,S2]. 



We now compare the behaviour with a type II vacuum solution. By (2.4) and 
( |14.7| ), we have 

(14.9) S'+ = -2{-^Nl + h,){Y.+ + 1) - ^(2 - 7)r!I]+ + ^TVf - 

-\ni{N2 + A^3) = \nI{2 - S+) + h^n + /i3, 

where \h^\ < 17e and |/i2| < 2 in [si,S2]. Let = {6 — 37)/4. Then, 



E+(s2) - S+(si) > a^{s2 - si)+ / 



{h2fl + hs.)dt. 



However, 



for all s e [si, S2], see (2.1). Thus, 

I / h2^lds\ < -n{s2)e 

We get 



S+(s2) - E+(si) > a^{s2 - si) - iee4(^^-^i) - 17e(52 - si). 

This inequality contradicts the statement that S2 — si may be taken equal to 4/a^, 
by choosing e small enough. We conclude that §2 — si < 4/a^ = ci.^, and that we 
may choose 02,7 = exp(16/a^). □ 

The following lemma deals with the decay in that has to follow an increase. The 
idea is that if Nx is on the boundary between big and small, and its derivative is 
non-negative at a point, then it will decrease as we go backward, and the solution 
will not move far from the Kasner circle until one of the other Ni has become large. 
That takes a long time and 17 will decay. 

Lemma 14.4. Consider a Bianchi IX solution such that 2/3 < 7 < 2. There is 
an e > such that if 

(14.10) iViiV2 + N2N3 + N3N1 < e 

in [si,S2], 

^iV2(.2) = ^(6-37), iV{(.2)>0 

and ri(s2) < C2^'ye, where C2.'y is the constant appearing in Lemma 14.!\ , then Vl 
decays as we go backward starting at S2, until s — si, or we reach a point s at 
which 

"(s) < 15 • 
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Proof. We begin by assuming that e > is a fixed number. As the proof progresses, 
we will restrict it to be smaller than a certain constant depending on 7. We could 
spell it out here, but prefer to add restrictions successively. Let A^i > e-^/^ in [^1,52] 
and iVi(ii) = e^/'* or ti = si, in case A^i does not attain e^/* in [si,S2]. As in the 
proof of Lemma 14. E , we conclude that Ni < e^^^, i = 2, 3 in [^i, S2], and that we 
may assume 



(14.11) 



si < 1. 



The variables (fi, S+, S_) have to belong to the interior of a paraboloid for N[ to 
be negative. Since N[{s2) > we are on the boundary or outside the paraboloid. 
The boundary is given by 5 = 0, where 

1 



9 



-(37 - 2)n + 2S^ + 2Ei - 4E+. 



An outward pointing normal is given by Vg, where the derivatives are taken in the 
order: and E_. Let 

£ = {Te [h, S2]:te [t, S2] =^ N[{t) > 0, nit) < C2,-ye}. 



Let T £ £. By ( 14.11 ) we get ^(t) < 2 and, as we are also outside the interior of 
the paraboloid, S+(t) < 1/2. For e, and thereby fi, small enough depending only 
on 7, we have 

cf. (|ll9|) . Using the above observations, we estimate in r, 



v<?-(r!',i]V,i]'_ 



.1/2 



where C-y only depends on 7. For e small enough, the scalar product is negative. 
Thus, if {^{t), I]+(r), E„ (t)) is on the surface of the paraboloid, the solution moves 
away from it as we go backward, so that iV{ > in [s,t] for some s < t. If we 
are already outside the paraboloid, the existence of such an s is guaranteed by less 



complicated arguments. As in the proof of Lemma 14.2, we get Q' > for e small 
enough depending only on 7, so that £ is open, closed and non-empty. Thus A^i 
decreases from S2 to ti going backward. Now, 



Si > 1 



1 



n-hi> -(37 + 2) - C2,^e - 3e 



in [ii, S2], so that 
(14.12) 



by an argument similar to Lemma 14.1, if e is small enough. We can assume e is 



small enough that the time required for A^i to decrease toe^ is great enough that 
if ti 7^ si, then the conclusion of the lemma follows by ( 14.12 ). □ 

Proof of Proposition Assume e is small enough that all the conditions of 



Lemma 14.2- 14.4 are fulfilled. We divide the interval [ti,T2] into suitable subinter- 
vals, such that we may apply the above lemmas to them. If 



(14.13) 



-iV2< 1(6 
4 ' - 8^ 



37) 



in T2 for i = 1, 2, 3, then we let ^2 G [''"ii '''2] be the smallest member of the interval 
such that (14.13) holds in all of \t2, T2]. Otherwise, we chose ^2 = T2- Either t2 = n 
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or 3iV^(t2)/4 > (6 — 37)/8, by a suitable permutation of the variables. If t2 7^ ti, 
let ti be the smallest member of [Ti,f2] such that 3iV^/4 > (6 — 37) /8 in [ti,i2]- 



Because of Lemma 14.2| , il d ecays in [^2,12]. If t2 — ti, we are done; let = 1. 



Otherwise, we apply Lemma 14.3 to the interval [ti,t2] to conclude that ^{t) < 
02,7^2(72 )in [ii, T2]. If ti = Ti, we can choose — 02,7. Otherwise, we apply 
14.4 to [ti,<i]. Either fl decays until we have reached ti, or there is a 



Lemma 



point si e [Ti,ii] such that fl{si) < VI{t2)/2. By the proof of Lemma 14.4 



assume that T2 — si > 1; some time has to elapse for the decay to take place. 

Given an interval [ti,T2] as in the statement of the proposition, there are thus two 
possibilities. Either ^{t) < C2.'y^l{T2) for all r € [ti,T2] or we can construct an 
si e [ti,'''2] such that T2 — si > 1, fl{si) < n(T2)/2, and fl{T) < C2.j^{t2) for all 
re [si, T2]. If the second possibility is the one that occurs, we can apply the same 
argument to [ti,si], and by repeated application, the proposition follows. □ 

Corollary 14.1. Consider a Bianchi IX solution with 2/3<7<2. // 

lim {N1N2 + N2N3 + N1N3) = 

r — ^ — 00 

and there is a sequence — > — 00 such that fl{Tk) — s- 0, then 

lim n{T) = 0. 

T — ' — 00 



15. Generic attractor for Bianchi IX solutions 

In this section, we prove that for a generic Bianchi IX solution, the closure of the 
set of type II vacuum points is an attractor, assuming 2/3 < 7 < 2. What we need 
to prove is that 

lim {n + N1N2 + N2N3 + N1N3) = 0, 

r — > — 00 

since then we may for each e > choose a T such that at least two of the Ni and 
must be less than e for t < T. The starting point is the existence of a limit point 



on the Kasner circle for a generic solution, given by Proposition 13.1 . Since there 
is such a limit point, there is a sequence Tk ~f ~oo such that Ni(Tk) and Q{Tk) go 
to zero. If 

(15.1) h = N1N2 + N2N3 + N1N3 

does not converge to zero, it must thus grow from an arbitrarily small value up 



to some e. By choosing e so that Proposition 14.1 is applicable, we have control 
over D,. A few arguments yield the conclusion that we may assume that it is 
the product A^2-^3 that grows, and that the growth occurs close to the special 
point (E+,E_) = (-1,0). Close to this point, f2, A^i and A^i(A^2 + N3) decay 
exponentially, so as far as intuition goes, we may equate them with zero. We 
thus have a Bianchi VIIq vacuum solution close to the special point (—1,0). The 
behaviour of N2N3 will be oscillatory, and we may reduce the problem to one in 
which the product behaves essentially as a sine wave. However, by doing some 
technical estimates, one may see that one goes down going from top to top during 
the oscillation, and that that contradicts the assumed growth. Figure |^ illustrates 
the behaviour. It is a simulation of part of a Bianchi VIIq vacuum solution. 
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Figure 7. Part of a Bianchi VIIq vacuum solution. 



We begin by rewriting the solutions in a form that makes the oscillatory behaviour 
apparent. Consider a non Taub-NUT Bianchi IX solution in an interval such that 
— 1 < S+ < 1. Define the functions 



(15.2) 



(15.3) 



y 



(1- 

V3 



N2 



1/2 



2 (1 



I]^)V2- 



The reason why these expressions are natural to consider is that, for reasons men- 
tioned above, TVi, and so forth may be considered to be zero. In the situation 
we will need to consider N2 — and E_ will have much greater derivatives than 
S+, so that it is natural to consider x and y as sine and cosine, since the constraint 
essentially says x'^ + y'^ = 1. Let 

(15.4) g - -3(A^2 + Ns) - 2(1 + S+)£y - gi + 52- 

In our applications, gi will essentially be constant, and g2 will essentially be zero. 
Lemma 15.1. The vector {x,yY satisfies 

x' = Ax + e, 

where A is defined as in (12.1), with g as in (15.4) e — {ex,ey)^, where the 
components are given by (15. i) and (15. t). 
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The error terms are 

(15.5) = 3N,y + (^7Vi(7Vi - N2 - N3) - ^{2 - - 
and 

(15.6) ey = [1(37 - 2)f^(l + S+) + ^(2 - 7)^^ + ^A^i(^i ~ N, ~ ^s)] 

It is clear that if we have a vacuum type VIIq solution, = = 0, so that we 



may write x = (sin(^(T)), cos(^(r))), where ^ is as in ( 12.2 ). In our situation, there 
is an error term, but by the exponential decay mentioned above, it only makes the 
technical details somewhat longer. 

We begin by proving that we can assume that the growth occurs in the product 
N2N3, and that fl can be assumed to be negligible during the growth. We also put 
bounds on S]-|_. They constitute a starting point for further restrictions. The values 
of certain constants have been chosen for future convenience. 

The lemma below is formulated to handle more general situations than the one 
above. One reason being the desire to prove uniform convergence to the attractor. 



We will use the terminology that if x constitutes initial data for ( p.l[ )-(2.3), then 
S+ (r, x) and so on will denote the solution of the equations with initial value x 
evaluated at t, assuming that r belongs to the existence interval. We will use 
<I>(t, x) to summarize all the variables. The goal of this section is to prove that the 
conditions of the lemma below are never met. 

Lemma 15.2. Let 2/3 < 7 < 2. Consider a sequence xi of Bianchi IX initial data 
with all Ni > and two sequences s; < r; of real numbers, belonging to the existence 
interval corresponding to xi, such that 

(15.7) lim d{Ti,xi) = 0, 

where d = n + N1N2 + N2N3, + N1N3, and 

(15.8) hisi,xi)>S 

for some S > independent of I. Then there is an e > and a ko, such that for each 
k > ko there is an Ik, a symmetry operation on $(-,a;;^), and an interval [uk,Vk] 
belonging to the existence interval of ^{-jXi^), such that the transformed variables 
satisfy 

{N2Ns){uk,xiJ = e, {N2N3){vk,xi,) < ee-^o^ eg-^o'^-i < {N2N3){t,xi^) < e 

(15.9) iVi(T,XiJ < eexp(-30fc) and 2 > A^2(r,a;iJ, Ns^t^xiJ > eexp(-25A:) 
for T ^ Furthermore 

(15.10) xi^) < e-"'^ and - 1< S+(-, J < 
in [uk,Vk]. 
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Remark. Observe that for the main apphcation of this lemma, the sequence xi will 
be independent of I. 

Proof. By (15.7) and ( 15. 8| ), there is an e > such that for every k there is a 
suitable and < Vk with [uk,Vk] Q [siki'''ik] such that 

(15.11) e-^°''-h<h{T,xi,) <2e 

h{uk, xif, ) = 2e, h{vk, xi^. ) = exp(— 20A:— l)e where r G [uk, Vk\- We can also assume 
that 

(15.12) h{T,xi,)<2e 
for all T G [uk^Tif.]. Furthermore, we can assume 

(15.13) (7Vi7V2A^3)(-,2;jJ <e2exp(-50fc-l)/4 

in [zifcjT/^]. The reason is that d{Ti,xi) converges to zero, so that (A^iA^2^3)(t;, a;;) 
also converges to zero. Consequently, we can assume {NiN2N3){ti^.,xi^.) to be as 
small as we wish, and thus we get ( |15.13| ) by the monotonicity of the product. 
Since we may assume r2(T;^,a;;^) to be arbitrarily small by (15.7), we may apply 
Proposition 14.1 in [uk,Ti^] by (15.12), choosing e small enough. Thus we may 
assume fl < exp(— 13fc) in [uk,Vk]. From now on, we consider the solution $(-,a;;^) 
in the interval [uk,Tii,] and only use the observations above. To avoid cumbersome 
notation, we will omit reference to the evaluation at xi^. By (15.11) and (15.13), 
we have in [uk, Vk] 



-20fe-l 



< ft = iViiV2iV3(— + — 



so that 



1 



1 

iV2 



1 



1_ > 4 



±) 



At a given r G [uk,Vk], one Ni, say A^i, must be smaller than eexp(— BO/c). If the 
second smallest is smaller than eexp(— 25fc), the largest cannot be bigger than 2, 
by Lemma ^.3| , but that will contradict h > eexp(— 20A: — 1) if fc is great enough. 
Thus, if A^i is the smallest Ni for one r, it is always the smallest. We may thus 
assume 

Ni < eexp(-30fc) and N2, N3 > eexp(-25fc) 



in [ufc, Vk]. If e is small enough, we can assume N2, -/V3 < 2 by Lemma 3.3. Thus, 
e-20fc-ig _ 4ee-^0k < j^^j^^ <2e + 4ee-30fc. 

We may shift Uk by adding a positive number to it so that 

(15.14) iN2m)iuk) = e and (7V2A^3)(t) < e 

for r G [uk,Vk]. We may also shift Vk in the negative direction to achieve 

iN2N3){vk) < ee-20fe, (iV2Af3)'(«fc) < and (A^2iV3)(r) > ee-^ofc-i 

for r G [uk,Vk]. The condition on the derivative is there to get control on S+. 

We now esta blish rough control of E+. Since (A^2^^3)'(wfe) < 0, -1 < T,+ {vk) < 0. 
Due to ( |15.9| ), ( [2!^ ) and the constraint, E'^ < if E+ = or S+ = -1. In other 
words, I]+(w/c) ~ implies > in [uk^Wk]. But if Uk < Wk then I]+(Mfe) > 
so that {N2N3){uk) < (^2-^3)(u'fc), contradicting the construction as stated in 
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Table 2. Subdivision of the interval of growth. 



Interval Bound on r 

[vk,(yk\ -4fc<r<-2fc 

[Tfe,Tfe] -4fc<r<-2fc 
[Tfc,rfc] -4fc < r 



(15.14). We thus have S+ < in [uk,Vk]- We also have — 1 < in that interval. 



□ 

Below, we will omit reference to the evaluation at xi^. to avoid cumbersome notation, 
but it should be remembered that we in general have a different solution for each 
k. Let 

r{T)= J^\q/2 + i:+)ds. 

Here we mean q{s,xi^) when we write q, and similarly for Observe that r 

depends on k, but that we omit reference to this dependence. All the information 



concerning the growth of N2N3 is contained in r, see (2^), and this integral will be 



our main object of study rather than the product N2N3. Let [ufc, Vk] be an interval 



as in Lemma 15.2 . Since 

(iV27V3)K) = e4'-("'=)(A^2A^3)K), 

we have r{uk) < — 5/c. Let Uk < i^k l£ '^k l£ Tk < Vk < Vk- Starting at Uk, let Vk 
be the last point r ~ — 4/c, so that r > — 4fc in [i^k, Wfe]. Furthermore, let r > —k in 
[fk, Vk] and finally, assume r < ~2k in [i^k, Tk]- We also assume that r evaluated at 
Tk, Tk, (Tk and Vk is —k, —2k, —3k and —Ak respectively. See Table [l5|. Why? The 
interval we will work with in the end is [ak, Tk], but the other intervals are used to 
get control of the variables there. First of all, we want to get control of and 
the interval [uk, Vk] together with the additional demand on Vk serves that purpose. 
The intervals at the other end, together with the associated demands, are there to 
yield us a quantitative statement of the intuitive idea that CI and A^i are negligible 
relative to the other expressions of interest. Finally, we need to get quantitative 
bounds relating the different variables; as was mentioned earlier, the main idea is to 
prove that N2N3 oscillates, but that it decreases during a period. In order to prove 
the decrease, we need to have control over the relative sizes of different expressions, 
and [vk, Cfc] is used to achieve the desired estimates. 



From this point until the statement of Theorem 15.1 , we will assume that the condi 



tions of Lemma 15.2 are fulfilled. We will use the consequences of this assumption, 
as stated above, freely. 

We improve the control of S+. Let us first give an intuitive argument. Observe 
that under the present circumstances, the solution is approximated by a Bianchi 



VIIq vacuum solution. For such a solution, the function Z-i, defined in (13.2), 
is monotone increasing going backwards. According to the Bianchi VIIq vacuum 
constraint, Z_i is proportional to (1 — E^)/Af2A^3- However, we know that N2N3 
has to increase by a factor of 6^°*^ going from Vk to Uk, and consequently 1 — 
has to increase by an even larger factor. The only way this can occur, is if a large 
part of the growth in N2N3 occurs when E+ is very close to —1. Taking this into 
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account, we see that the relevant variation in 1 — = (1 — S+)(l + S+) occurs in 
the factor 1 + Below, we will use the function (1 + Yj+)/N2N^ instead of 
Let us begin by considering the vacuum case, in order to see the idea behind the 
argument, without the technical difficulties associated with the non- vacuum case. 
We have 



(15.15) 



1 



< 



in our situation, cf. Lemma 15. 3| and ( 15.10 ). For r G [i^k, Vk] we get 

(iV2^3)(r) 



0<l + I]+(r)<(l + S+(ufe))- 



< e 



-4k 



iN2N3)iuk) 

by our construction. 

Let us make some observations before we turn to the non-vacuum case. First we 
analyz e the derivative of (1 + I]_|_)/A^2-^3 in general. The estimates ( 15.16| ) and 
( 15.17| ) will in fact be important throughout this section. 



Lemma 15.3. 

(15.16) 

and 
(15.17) 



Let Uk and Vk be as above. Then 
■l + S+\' , -2[(1 + S+)2 + S^](l + !]+) + 1(2 -7)f7 



N2N3 



< 



in the interval [uk,Vk] for k large enough. 



Remark. Observe that 1 + S_|_ > in [uk,Vk] by (15. IC), so that the first term 
appearing in the numerator of the right hand side of (15.16) has the right sign. 

Proof. Using (^.4|), we have 



1 



= [-(2 -2n- 2E^ - 2S^)(S^ 



l)--(2-7)r!I]++ 



+ -N,{N, 



N2 - N3) - {2q + 4S+)(1 + j:+)]{N2N3y 



Consider the numerator of the right hand side. The term involving the Ni has the 
right sign by ( 15.91 ), and the terms not involving add up to the first term of the 
numerator of the right hand side of ( 15.161 ). Let us consider the terms involving Q. 
They are 



2n{i + s+) - ^(2 - 7)17(1 + E+) + ^(2 - j)n - (37 - 2)n{i 

= -i(37 - 2)r!(l + E+) + ^(2 - 7)f7 < ^(2 - j)n 



s+) 



proving (15.16). To prove (15.17), we observe that by the constraint and the fact 



that < 1 + I]+ < 1 in the interval of interest, we have 

-(2 -2n- 2S^ - 2E^)(E+ + 1) < 3Ni{N2 + N3){1 + E+) < 37Vi(iV2 + N3). 



Inserting this inequality into (2.4), we get 



S'+ < -^(2 - iMl + S+) + ^(2 - 7)r! + ^N,i9N, - 3iV2 - STVg) < ^(2 - l)^ 
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by (15.9) and (15. IC) if k is large enough, proving (15.17). □ 



In the vacuum case, is monotone in our situation, see (15.17), but in the general 
case we have the following weaker result. 

Lemma 15.4. Consider an interval [s,t\ C [uk^Vk] such that 



Then 
(15.18) 

if k is large enough. 
Proof. In [s,t] we have 



St> g(37 + 2) 
(l + E+(<))-f^(<) < 1 

f2' > a^rj. 



E+(s) 



where — 3(2 — 7)/2, see the proof of Lemma 14.1 . Thus, 

n{u) < n{t) e-Kp[a^{u - t)] 
for all u e [s,t]. Integrating ( |l5.17| ) we get ( |15.1S| ). □ 
In connection with (15. 1£), the following lemma is of interest. 

Lemma 15.5. If k is large enough and 



for some r G [ufc,ffe]; then 



(1^ 
(1 



S+)3 



> 4(2-7)^^ 



m [Uk,T\ 



Proof. If the solution is of vacuum type the lemma follows, so assume > 0. Let us 
first prove that (l + E+(it))^ > e'^f2(r) for u G [uk, t]. Assume there is an s e [itfe, t] 
such that the reverse inequality holds. Then there is a t with t > t > s, such that 
(H-E+)^ < e^''fl(r) in [s,t], with equality at t. Because of ( |15.10| ), Lemma p^^i] is 
applicable for k large enough. Thus 

(15.19) e'^n^/^ir) - n{t) < 1 + E+(s) < e'^^^n^/^ir). 

However, by the proof of Lemma |15.2| , Proposition 14.1 is applicable in any subin- 
terval of [uk,Vk], so that n{t) < Cjil{T). Substituting this into (15.19), we get 

e'^T!i/3(r)-c^r!(T) <e'=/3^i/3(r), 
which is impossible for k large enough. 
Thus we have, for u € [uk, t] and k large enough, 



(l + E+(w))3 > e'=17(r) > e* 



1 



(2- 7)2^4 4 



;/^7 

where is the constant appearing in the statement of Proposition 
follows. □ 

We now prove that wc have control over 1 + in [vk, ^^fc]- 



14.1 



The lemma 
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Lemma 15.6. Let Vk o,nd Vk be as above. Then for k large enough, 



(15.20) 

in [vk^Vk]. 



< 1 



< e 



Proof. Assume 1 + S+(t) > e for some t € Because of ( 15.1C| ), we then 

conclude that Lemma 15.5 is appHcable, so that 



1 



< 



in [ufe,r] by ( |l5.16| ). Thus 



> 



1 + S+(t) 



> 



{N2m)iuk) - {N2N3){t) - {N2N3){t)' 
but by our construction 

(A^2A^3)(t) = e4'^("'=)-4'-(-)(7V27V3)(^^fe) < e-20'=+i6fe(^r2Ar3)(y,), 

so that 



e^*^ < 1 



S+(Mfe) < 1. 



The lemma follows. □ 

Corollary 15.1. Let and Vk be as above. For k large enough, 

n + + (l + E+)2 <4:e-'' 

in [vk,Vk]- 



Proof. By (15.£), we have 



Ni{N2 +N3) < Aee 



-30k 



in [uk,Vk]- This observation, the constraint, and Lemma 15.6 yield 



in [vkTVk], for k large enough. The corollary follows using Lemma 15.6 



□ 



The 
S+ is 



The next thing to prove is that Ni and Q are small compared with 1 + S 
fact that r{rk) — —k will imply that the integral of 1 + S+ is large, but if 1 
comparable with A^i or 17, it cannot be large since Ni and decay exponentially. 

The reason (1 + appears in the estimate ( 15.21 ) below is that the final argu- 
ment will consist of an estimate of an integral up to 'order of magnitude'. Expres- 
sions of the form {1 + and {1 + {N2 + will will define what is 'big' 
and 'small', and here we see to it that terms involving and Ni are negligible in 
this order of magnitude calculus. Finally, the factor exp(— 3fc) is there in order for 
us to be able to ignore possible factors multiplying expressions involving Ni and 
51. We only turn up the number k and change exp(— 3fc) to exp(— 2fc) to eliminate 
constants we do not want to think about; consider (15.5) and (15.6). 

Lemma 15.7. Let Vk o-iT-d- Tk be as above. Then for k large enough, 

(15.21) n + Ni + Ni{N2 + N3) < e"3'=e3''-'(^~"'='(l + i;+)^ 
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in [t'fcjTfe] where b-y > 0. Furthermore, 



(15.22) 
in [wfc,Tfc]. 
Proof. Note that 
(1 

so that 
(15.23) 

Let 



1 



N2N3 



(1 



N2N3 



nVk rVk rVk 

/ (l + S+)dr< / (E^ + S+)dr< / (g/2 + S+)dT 

J n. J Tk J ri- 



k< I (1 + S+)dr. 



pi^r! + iVi+7Vi(7V2+iV3). 



By the construction in Lemma 15.2 



we may assume 

12fc 



Pi{vk) < e 



Because of Corollary 15.1, we have 

Pi{t) < e-i2fee4^(--"'=) 
for all T G [i>k, Vk], where 6^ > is some constant depending only on 7. Let 
P2(t) = e-^'^e^'^^^-"") > e^'^e-^''^ Vi(^)- 



The assumption that (H-S+)^ < p2 in [r^, u^] contradicts ( 15.23 ). Thus there must 
be a to G [?'fc,i'fc] such that (1 + S+(io))^ > P2{to)- In the vacuum case, 1 + E+ 
increases as we go backward, and p2 obviously decreases, and thus we are in that 
case able to conclude (1 + S+)^ > p2 in [j^fc,r/j]. In the general case, we observe 
that (1 + £+(to))^ > e.^^^i'to) by the above constructions. We get 



1 



N2N3 



< -2E! 



(1 + S+) 
N2N. 



in [ufe,to]7 by combining Lemma 15.5 and ( 15.16 ). Inequality ( |15.22| ) follows. Thus, 
if T e [i/fe, Tfe], we have 

1 + I]+(r) > + S+(to)) > e^^il + S+(to)). 

Consequently, we will have (1 + I]+(r))^ > P2{t), since 1 + has increased from 
its value at smd p2 has decreased. The lemma follows. □ 

Next we establish a relation between 1 + E+ and the product N2N3. We prove 
that (1 + S+)/(A^2^3) can be chosen arbitrarily small in the interval [ak,Tk], by 
estimating it in Vk, and then comparing the integral of 1 + S+ from to ak with 
the integral of E'i over the same interval. The following lemma is the starting 
point. 

Lemma 15.8. Let (Tk, Tk be as above. Then for k large enough, 



(15.24) 



1 + E+(t) 1 , „ 
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z/tG [tTfcjTfe]. Furthermore, 



l + S+(r) ^ 1 
(A^2A^3)(t) - e 



Proof. The statement follows from ( |15.22|) , and the fact that 



(i + s+)K) 

{N2N3){uk) 



1 

< -. 



□ 



Considering the constraint, it is clear that should be comparable with 1 + 
when N2 — N 3 and E_ oscillate, and thus the integral should be comparable with 
fc, cf. ( |l5.23|) . However, we have to work out the technical details. 

We carry out the comparison between the i ntegr als in three steps. First, we estimate 
the error committed in viewing x and y in ( 15. 2| ) and ( 15.3 ) as sine and cosine. Then 
we may, up to a small error, express the integral of S?. as the integral of sin^(?7/2), 
multiplied by some function /(??) by changing variables. In order to make the 
comparison, we need to estimate the variation of / during a period: the second 
step. The only expressions involved are 1 + E+ and N2 + N3. The third step 
consists of making the comparison, using the information obtained in the earlier 
steps. 

Let X, y, g, gi and (72 be defined as in ( 15.2[ )-( p^!4| ), and ^, x and y be defined as in 
the statement of Lemma 12.1, with tq replaced by and by (pk- Observe that 



a;, y and ^ in fact depend on k. We need to compare x with x. 
Lemma 15.9. Let Vk cind be as above. Then for k large enough, 



(15.25) \Y? 

in[vk,Tk\. Furthermore, 

(15.26) 

and 

(15.27) 

in that interval. 



< 12e-2fe(i + 



■f) 



< e~ 



|x-x|| < 3e-2fe(i 



Proof. We have 

II - mr.) . fir.)m . II - fl + i^^(^^^^3)-!A^f-» 



(15.28) 



< e 



-2fc 



(l + S+(r,))« 



1/2 



< 



by (|15.21| ). Equation ( |l5.26| ) follows similarly. By (|l5.5| ), ( |15.6D , ( [l5.2l[ ) and 
( 15.26| ), we have 

||e(.s)|| < 2^6-2^(1 + S+(.s))«e3''-(*-"^) 

for k large enough. Let us estimate how much 1 + S4. may decrease as we go 
backward in time. By (15.17) and (15.21), we have 



(1 + S+)' < |(2 - 7)e-3fee3''T(--"^-)(i + £+)« 
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SO that if [s,t] C [z^fc,Tfe], 

(15.29) 1 + E+(t) < exp(exp(~2/c))(l + S+(s)), 
for k large enough. Thus, for t < t/c, we get 

(15.30) l|e(s)||ds < e-2'=(l + S+(t))8. 

By (|l2j), ( |15.30D , ( |15.29| ) and ( |15.28| ), we thus have 

5 



!x-x|| < -e 



-2k/ 



'•(i + s+)« 

in [i^k,Tk], and ( 15.27 ) foUows. Since < 1 and \x\ < 1.1, of. ( 15.26| ), we have 

-x^\<6e-^''il + ^+f, 

so that 

- (1 - S^)a;2| < 12e-2fc(i + 

in the interval [j/fe,rfc]. □ 
Let us introduce 



(15.31) 



77(t) = 2^(t) = 2 / 

Jti, 



g{s)ds + 2(t>k., 



where g = —3{N2 + N^) — 2(1 + S+)a;?/ = gi + g2- The reason we study rj instead 
of ^ is that the trigonometric expression we will be interested in is sin^(i^), which 
has a period of length tt, cf. Lemnia [l5.9| . In the proof of Lemma 15.10, it is shown 
that, in the interval [i^fe,Tfe], the first term appearing in g is much greater than the 
second. We can thus consider functions of r in the interval [vk, Tk] to be functions of 
77. We will mainly be interested in considering an interval [770, ?7o + 27r] at a time, so 
that we will only need to estimate the variation of the relevant expressions during 
one such period. 

Lemma 15.10. Let rji^k = vif^k) and 772, fc = vi'^k)- U + 27r] C [?7i,fc, ?72,fc] 

and ria,rii) G [771, 771 + 27r], then for k large enough 



(15.32) 

(15.33) 

and 
(15.34) 



-6ir/e ^ 



1 ^ l + ^+iVa 



< e^'"l\ 



2 - l + S+(77b) 



< 2 



Proof. Because of Lemma 15.8 , 

1 + S+ l + S 

(15.35) — ^ < 



N2 + N3 - 2{N2N3y/^ 
1/2 



l5i|/2< IffI <2|5i|. 



< 



1 



N2N3 



2e \ iN2N3){uk) 



{N2N3)^/^{uk) < 



2N2N3 
1 



< 



-2k 



2ei/2 
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in the interval [vk , 



By (15.26) we may assume x'^+y'^ < 2 in [i^k, Tk]- Combining 
this fact with (15.35) yields (15.34) in [j/fe,rfe]. Thus, drj/dr < in that interval. 
We have 



d{N2 + N3) 



\^iiq + 2S+)(A^2 + Ns) + 2V3I]_(A^2 - A^s))! < 



1 



< ^(37 - 2)n + + (1 - + S+l + 2^(1 - S^) < 6(1 + S+ 

2 |.g| 

so that 

, 1 d{N2 + N3) 



, 1 + 



N2 + N3 dri 



< 6 



1 + S4 



1 + 



< 6 



+ 2- 



iV2 + TVs (iV2 + iV3)2 - 7V2 + 7V3 A^2A^3 



< - 



in [vi;, cTfc] for k large, by Lemma 15.8 and (15.35). If N2 + has a maximum in 
?7max e [vijVi + 27r] and a minimum in r^min, we get 

(iV2+iV3)(77max) ^ err/e 
{N2 + m){7Jnun) - 

and (|15.32| ) follows. We also need to know how much 1 + S+ varies over one period. 
By is 

(1 + = (2E^ + 2I]2_ - 2)(1 + E+) + /i, 

where /i is an expression that can be estimated as in ( 15.2l| ), so that we in [i^k, Tk] 
have 



' l + E 
for k large enough. Thus, 

(15.36) 



< 2(1 - S;)(l + x') + (1 + S+) < 13(1 + 



1 d(l + ^ 10(1 + 1]+) 



' 1 + S+ d?7 ' - N2 + N3 
so that ( 15.33| ) holds if k is big enough and \Va - Vb\ < 27r by ( |15.35| ). □ 
Lemma 15.11. Let (ifc anc? be as above. Then if k is large enough, 



N2N3 - e 



in [(7k, Tk] where q > 0. 



Proof. Observe that similarly to the proof of Lemma 15.7, we have 

(1 + S+) 



k< / {l + Y.+ )dT 



-2g 



-dr]. 



The contribution from one period in r/ is negligible, by (15.35) and ( 15.34| ). Compare 
this integral with 



'rii,k ~5 

Now, 



V2.k (1 _ 5]2_)^2 



drj 



VI, k 



V2,k _ (1 - Ei)a;2 



ni.k 



-9 



-drj ^ Ii^k + h,k- 



hA < e- 



l + S^ 



■dri 
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by (15.25). Consider an interval [?/i,?7i + 27r]. Estimate, letting rja and rjh be the 
minimum and maximum of respectively, and r/min, Jyinax the min and max for 
gi in this interval, 



drj > 



V1+27T 



(l + S]+)x^ 

-5 



drj = 



?ji+27r 



(l + S+)sin2(,7/2) 



-9 



drj > 



2|.gi(?7inax)| ~ 2 Iffll'ymin)! ~ 4 |gi(?7min)| 



1 + 



1)1 



>)l+27r 



l + S+(7?) 



d?7. 



IsiC??.™.!)! 

where we have used ( 15.32[ ), ( 15.33| ) and (15.34). Assuming, without loss of gener- 
ality, that 772, fc — ?7i,fc is an integer multiple of 2tt, we get 

2Y?_dT = / d-q = h,u + /2.fe > — 



l + S+(77) 



1 

> — e 
- 20 



— 67r/e 



l + S+(r7) 
-5'('7) 



1 

10' 



/ (1 + S+)dT > = cfc 



dry > 



for fc large enough and the lemma follows from (15.24). □ 

The following corollary summarizes the estimates that make the order of magnitude 
calculus well defined. 

Corollary 15.2. Let ak and Tk be as above. Then 

1 + <i,-.. 



(15.37) 

(15.38) 

and 
(15.39) 



{N2 + A^3)' 
1 + S+ 



N2 + NS, 



< e 



-2k 



1 _ e-2fc < A < 1 + g-sfe 



91 



in [crfe,Tfc] for k large enough. 



Proof. Observe that by Lemma 15.11 



1 + ^+ < < ig-c.. 

{N2 + N^)^ - N2N3 - e ' 



and 



1 + E+ ^ 1 + S+ 



1 



iV2 + iV3 - 2(iV27V3)i/2 - ' 2e -' 



2fe 



for fc large enough, cf. (15.35). We have 



9_ 

91 



2(l + S+)iy 



3(iV2 + N3) 

By (15.26) and the above estimates, we get ( 15.38| ) for k large enough. □ 

The interval we will work with from now on is [cTfc, t^]. Let 77 be defined as in (15.31), 
but define r]i_k — vi'^k) and 772, ^ = rj{ak). We need to improve the estimates of the 
variation of 1 + and N2 + -/V3 during a period contained in [771,^, 772,fc]- 
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Lemma 15.12. Consider an interval 2 ~ + 2tt] C [ryi^fc, ?y2,fc]; where rji^k = 

ri{Tk) and r]2.k = vi<^k)- Let rja and r/h correspond to the max and rain of 1 + in 
2, and let Ty^ax o,nd Tymin correspond to the max and min of N2 
interval. Then, 

407r(l + E+(77fc))2 



A^3 in the same 



(15.40) 
and 
(15.41) 



|E+(r;fa)-E+(r,,)| < 



(7V2+iV3)(?7,„ax) 



) /207r 

y < exp( exp(— Cefc)). 



{N2 + N3){n 

Proof. The derivation of ( 15.36| ) is still valid, so that 

I 1 (i(l + S+)| ^ 10(1 + S 



drj 



N2 + N3 



By ( 15.38 ) we conclude that (1 + S+(?7q))/(1 + S+(ryf,)) can be chosen to be arbi- 
trarily close to one by choosing k large enough. Now, 



1 



d{N2 + N3) 



1 



1 diN2+N3) 



N2 + N3 



drj 



dr 



1 



1 



N2 + N3 2g 
and consequently 



((g + 2I] + )(iV2 + iV3) 



N3 2g 
2\/3I]_(A^2-A^3)) = 



q + 2E+ 4(1 - i:l)xy 



2ff 



2{N2 + N3)g' 



1 



d{N2 +N3). 



10 
< — 
e 



'N2 + N3 drj 

Equation ( 15.41 ) follows, and the relative variation of N2 + 
can be chosen arbitrarily small. Finally, 

l + S+(r?a) 



A^3 during one period 



\i:+irjk)-^+iVa)\ = {l + ^+iVb))\ 
307r(l + 



< 



1 - 



s+(j?&) 



-II < 



(A^2 + iV3)(r;^i„) 

We may also change ri„ 



to 77inax at the 



by (15.36) and the above observations, 
cost of increasing the constant. □ 

As has been stated earlier, the goal of this section is to prove that the conditions 
of Lemma 15.2 are never met. We do this by deducing a contradiction from the 
consequences of that lemma. On the one hand, we have a rough picture of how the 
solution behaves in [(Tfc,rfc] by Lemma |15.9| , Lemma 15.12 and Corollary 15. S. On 
the other hand, we know that, since r{ak) — t{t]^) = — fc, 

(15.42) 



(4(37 



2)f7- 



S+)(ir = a J. 



VI, k 



-2g 



-dr]. 



We will use our knowledge of the behaviour of the solution in [(7^, Tk] to prove that 
( 15.42| ) is false. Observe th at r]i, k < '>l2,k, and that the contribution from one period 
is negligible, cf. Corollary 15.2. Also, ^ as /c ^ cxd so that we may ignore 
it. We will prove that for k great enough, the integral of (E^ + + E+)/(— 2f/) 
over a suitably chosen period is positive. From here on, we consider an interval 
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[771 , 771 + 27r] which, excepting intervals of length less than a period at each end of 
[ili.k,V2.k], we can assume to be of the form [— 7r/2, 37r/2]. There is however one 
thing that should be kept in mind; when translating the ry-variable by 2m7r the 
^-variable is translated by rmr. In other words, there is a sign involved, and in 
order to keep track of it we write out the details. By the above observations we 
have. 



y37r/2+2m2,fcJr j-i2 1 j-i2 1 j-i 

(15.43) -k = (3k+ / -df], 



Lemma 15.13. For each k there are integers mi^k and m2.k such that 

-7r/2+2mi,fc7r ^"^9 

where (3k as fc — > oo, and 

r]i,k < -7r/2 + 2mi,fc7r < rji^k + 27r, r]2,k - 2tt < 3tt/2 + 2m2,fc7r < 772, fe. 

Consider now an interval 

[-tt/2 + 2mTT, 3tt/2 + 2m7r] C [-tt/2 + 2mi,fc, 37r/2 + 2m2,fc7r], 
where m is an integer, and make the substitution 

f] — rj — 2TO7r, ^ = ^ — TOTT 

in that interval. Compute 

+ (1 - i:\)x^ + E+ = (1 + S+)(E+ + (1 - E+)i(l - cosr;)) = 

= (1 + S+)(i(l + S+) - \{l - S+) cos 77) = 1(1 + E+)((l + S+) - (1 - S+) COS77). 
This expression is the relevant part of the numerator of the integrand in the right 



hand side of (15.43). There is a drift term yielding a positive contribution to 
the integral, but the oscillatory term is arbitrarily much greater by Lemma |15.6 



The interval [— 7r/2, 37r/2] was not chosen at random. By considering the above 
expression, one concludes that the oscillatory term is negative in [— 7r/2,7r/2] and 
positive in [7r/2, 37r/2]. As far as obtaining a contradiction goes, the first interval 
is thus bad and the second good. In order to estimate the integral over a period, 
the natural thing to do is then to make a substitution in the interval [7r/2, 37r/2], 
so that it becomes an integral over the interval [— 7r/2, 7r/2]. It is then important to 
know how the different expressions vary with 77. We will prove a lemma saying that 
E+ roughly increases with 77, and it will turn out to be useful that is greater in 
the good part than in the bad. Let 

3^/2+2m^ Y.l+T?_+Y.+ _l f^-^/^ (1 + I]+(77 + 2777,7r))2 



(15.44) J = / ^ ^d77 - - / ^ -i^^- ^d77 

1 /"^^/^ (1-S2 (r7 + 277Z7r))cos77 

-077+ 



2j_,/2 -25(77 + 2m7r) 

3t/2 Y,'i{ri + 2mTT) - (1 - i;+(77 + 277?,7r))2a;2(i7 + 277i7r) 

, r dfj = Ji + J2 + Js- 

-^/2 -25(77 + 27777r) 

If we can prove that J is positive regardless of ttt. we are done, since J positive 



contradicts (15.43). The integral Ji is positive, and because the relative variation 
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of the integrand can be chosen arbitrarily smaU by choosing k large enough, Ji is 
of the order of magnitude 

^ ' N2 + N3 

If negative terms in J2 and J3 of the orders of magnitude 



(15.46) 
or 

occur, we may ignore them by (|l5.38| ) and (|15.37D . By (|l5.25D , J3 may be ignored. 



(iV2 + Ns)^ 



Observe that the largest integrand is the one appearing in J2- However, it oscillates. 



Considering ( 15.44 ), one can see that writing out arguments such as fj + 2rm: does 
not make things all that much clearer. For that reason, we introduce the following 
convention. 

Convention 15.1. By S+(77) and S+(— 77 + tt), we will mean T,+ {fj + 2rrm) and 
!]+(— 77 + TT + 2m7r) respectively, and similarly for all expressions in the variables of 
Wainwright and Hsu. However, trigonometric expressions should be read as stated. 
Thus cos(7y/2) means just that and not cos(7y/2 + rrm). 

Definition 15.1. Consider an integral expression 

1= / /(ry)dry. 

J-7r/2 

Then we say that I is less than or equal to zero up to order of magnitude, if 

r.3ir/2 



where g satisfies a bound 



-7r/2 



g < Ci + 62 



(A^2+iV3)2 '(iV2+A^3)'' 

for k large enough, where Ci and C2 are positive constants independent of k. We 
write / < 0. The definition of / > is similar. We also define the concept similarly 
if the interval of integration is different. 

We will use the same terminology more generally in inequalities between functions, 
if those inequalities, when inserted into the proper integrals, yield inequalities in 
the sense of the definition above. We will write w if the error is of negligible order 
of magnitude. 

Lemma 15.14. If J2 as defined above satisfies J2 ^ 0, then J is non-negative for 
k large enough. 

Proof. Under the assumptions of the lemma, we have 

1 (iJ:S+)2^^_ /•3-/2^^ (i + s+)3 (i + s+)3 

-tt/2 J--^/2 



THE BIANCHI IX ATTRACTOR 



55 



-7r/2 



-2.9 



-dr]. 



By Corollary 15.2, Lemma 15.12 and (15.25), we conclude that for k large enough, 
J is positive. □ 

The following lemma says that S+ almost increases with fj. 
Lemma 15.15. Let —tt/2 < fja < fji, < 37r/2. Then 

S+(77fc) - S+(r),) > -(1 + I]+(r)„,i„))8, 
where Tymin corresponds to the minimum o/ 1 + S+ in [— 7r/2, 37r/2]. 

Proof. We have 

3 
2' 



so that 



— + > - 2-7) — . 
dfj - 2^ '^2g 



Using ( |15.21| ), (|15.38|) and Lemma |l5.12| , we conclude that 



+ 1 



The lemma follows. □ 
Lemma 15.16. // 

satisfies / ^ 0, then J2 ^ 
Proof. Consider 

-J9 - 



1 = 



1 + S+ 



37r/2 
-7r/2 -5 



cos fidrj 



.3^/2 (1_5]2)C0S7) 

^ dfj 

-7T/2 



■^'^/^ (I]+(3V2)~S+)(1 + S+) 



+ (1-E+(3V2)) 



-7r/2 
-7r/2 



cos f)df}-\- 



■ cos f]df). 



The first integral is negligible by (15.40). The lemma follows. □ 
Lemma 15.17. // 

(l + S+(77))(gi(r))-gi(-r) + 7r)) 
/-V2 9{ri)g{-fl + -n:) 



h 



cos fjdf) 



satisfies Ii < 0, then Ji > 0. 
Proof. We have 

^3^/2 -|_ ^ 



/ = 



.3^/2 1 + , 1 + s+ , /■^''/^ 1 + s+ 

/ cos r\dr\ = I cos 77077 + / cos 77^77 

J—n/2 -5 J-Ti/2 -9 J-KI2 -5 
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Make the substitution x = + in the second integral; 



-7r/2 



l + ^+i-X + n) 



Thus, 



-.9(-X + 7r) 



cos(-x + 7r)(-dx) 



7r/2 



l + I]+(-X + 7r) 
-7^/2 -5(-X + 7r) 



cos(x)dx- 



/ = 



-7r/2 



COS rydr) 



/2 (1 + S+(-r) + ^))5(7y) - (1 + I]+(77))5(-7; + ^) 



7r/2 



cos 77^77. 



But 



(1 + I]+(-7? + 7r))5(77) < (1 + S+(7?))5(?7) 
by Lemma 15.15, so that 

r-/2 (l_+S+(77))(5(77)-.9H; + 7r)) 

-7r/2 



(15.48) 
Now, 



/ < 



aiiDgi-v + tt) 



cos fjdfj. 



giv) - gi-v + tt) = gi(?7) - 5i(-?7 + tt) + 52(7?) - .g2(-7? + tt), 

but since 2xy = sin fj and the error committed in replacing x with x and y with y 
is negligible by (15.27), we have 

g2iv) - g2i~V + tt) « -(1 + S+(77)) sin 77 + (1 + S+(-77 + tt)) sin(-77 + tt) = 

= (S+(-^ + ^)) - S+ (7/)) sin 77. 

The corresponding contribution to the integral may consequently be neglected; the 
error in the integral will be of type ( 15.47 ) by ( 15.40 ). Consequently, if 

(^E+(77))(5i(77)-5i(-77 + 7r)) 

-tt/2 



h 



g{fi)g{-f] + 'K) 



■ cos ?7c?77 



satisfies Ii < 0, then / ^ by (15.48), so that the lemma follows by Lemma 15.16 



□ 

Let 



hi{v) = giiv) - gii-fj + Tr). 
We estimate hi by estimating the derivative. We have /ii(7r/2) — 0. 

Lemma 15.18. Let hi be as above. In the interval [— 7r/2, 7r/2], we have 



(15. 



Proof. Compute 



dhi 
dfj 



> 3 



1-S|(77) , l-S^(-77 + ^) 



-giv) 



-gi-fj + tt) 



sm 77. 



But 



dhi,.. dgi,.. dgi 

-priv) = -7^(77) + -pri-V + tt) 
dr] df] drj 



^ = -|-((g + 2S+)(iV2 + iV3) + 2V3S_(iV2-7V3)) = 
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25 9 
Observe that x and y arc trigonometric expressions, and that 

2x{f) + 2rmT)y{fj + 2rm:) = 2sin(7)/2 + tott) cos(r)/2 + 77i7r) = sin 77. 

We have 

V3I]_(7V2 - A^a) « 2(1 - i:\)xy = (1 - E^) sin 77, 

so that 

^ « (i(37 - 2)r! + + + E+) - ^(^(37 - 2)^! + + T.I + 

3(1 - S^) sin 77 
9 

The middle term and all terms involving may be ignored. Estimate 

S^(77) + Y?_{f,) + S]+(77) + S^(-7) + ^) + E2_(-7y + tt) + S]+(-7y + ^) « 

« S^(7)) + (1 - E2^(77))(sin2(77/2 + tti^) - 1/2) + hi - £^(77)) + 



+S^(-77 + tt) + (1 - S^(-r7 + 7r))(cos2(77/2 + r7i7r) - 1/2) + -(1 - S^(-7y + 7r)) + 



+E+(-77 + tt) = 1(1 + S+(77))2 + 1(1 + E+(-7y + ^))2+ 

+ (1 - S^(7y))(sin2(,7/2) - 1/2) + (1 - + ^))(cos2(77/2) - 1/2). 

The first equality is a consequence of ( 15.25| ). Due to the fact that 77 G [— 7r/2,7r/2], 
we have cos^(77/2) — 1/2 > 0. Since —fj + vr > 77 and Yj^ increases with fj up to 
order of magnitude according to Lemma 15.15, we have 

l-S2^(-77 + 7r)> 1-^2^(77). 

Consequently, 

i(l + I]+(r7))2 + i(l + E+(-77 + vr))^ + (1 - E^(7y))(sin2(7y/2) - 1/2) + 

+ (1 - + 7r))(cos^(7)/2) - 1/2) > 1(1 + Y+{ft)f + 1(1 + E+(-7) + ir)f + 

+ (1 - E^(r)))(sin2(r)/2) - 1/2) + (1 - S^(r)))(cos2(77/2) - 1/2) > 0. 



In other words, we have ( 15.49 ). Here the importance of the fact that S4. is greater 
in the good part than in the bad becomes apparent. □ 

Lemma 15.19. Let Ii he defined as above. Then Ii < 0. 

Proof. Let Ty^ax and T^niin correspond to the max and min of —5 in the interval 
[— 7r/2, 37r/2], and let fja and Tyf, correspond to the max and min of S+, in the same 
interval. Observe that for fj € [— 7r/2, 37r/2], we have 

l-E^ (77a) >1-S^ (77) >1-S^ (??;.). 
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In order not to obtain too complicated expressions, let us introduce the following 
terminology: 

1-T.l{fj) l-i:l{i^a) 

ai — 6 — ^ < 6 7— — < 6 ^ = 02 and 

-5'('7max) -gW -giVmin) 

^ l + S+(7?fc) ^ l + ^ l + ^+iVa) ^ 

^ .g^('7max) ~ g{'n)9{-f] + it) ~ 5^(77min) ^' 

where t) G [— 7r/2, 37r/2]. Observe that 

(15.50) lim — = lim ^ = 1, 

k — *C30 0,2 k — *oo 02 

by Corollary 15.2 and Lemma 15.12. Consider the interval [0,7r/2]. By ( 15.49| ), we 
have 

(15.51) — — !- > fli sinry, 

drj 

so that 

hi{fj) = hi{TT/2) — / — — d?? < — fli COS77 

Jf, d-n 

in the interval [0, 7r/2]. Now consider the interval [— 7r/2,0]. We have 

dhi 
dfj 

Consequently, 

hi{fj) = hi{0) - [ dfj < -ai + a2{l - cos fj) 
Jf, dri 

in the interval [— 7r/2,0]. Estimate 

r-/' {i + j:+{mgiiv)-gi{-v + ^)) 
g{f])g{-f] + Ti) 



> 02 sin 77. 



cos fidfi ■ 



r'^ (i+s+(77))fei(77) /-/^ {i+^+m , 2^,,-^ 

/ — .-s . s cosridri< / —— -— — -(-oicos 77)^77 < 

^ g{v)g[-v + '^) Jo g\mg(.-ri + iT) 

2 ~,~ iraibi 



.7r/2 

< —aibi / cos^ 7)^77 

Jo 1 

We also estimate 

'•° (1 + E+ (77)) (51 (77) (-77 + 71)) 



V2 g{v)g{-v + -^) 



cos 77(^77 



(l + S+(77))/zi(77) ^ (l + S+(77)) 

V +\uj I ' COS rjdr] < -ai / cosr]dr]+ 



+02 



^/2 g{fi)g{-fi + T^) ^ J~7:/2 giv)gi-fi + TT) 

(1 -{- Xj (")) 

, . , ^ -(1 — cos 77) 00377^77 < —0161 / COS fjdfi+ 

-7r/2 g[V)g{-V + T^) J -7^/2 

f° ^ 

+O2&2 / (1 - cos 77) COS 77(^77 < -Oi&i + (1 - —)0262- 
J--K/2 4 
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Adding up, we conclude that 

h < -(1 + V4)aifei + (1 - TT/A)a2b2 = [-(1 + 7r/4)^ + (1 - 7r/4)]a2fe2 

0202 



which is negative for k large enough by ( 15.50 ). Thus /i < 0. □ 



Theorem 15.1. The conditions of Lemma 15. i are never met 



Proof. If the conditions are met, then Lemma 15.13 follows, and also that it is false 



by Lemmas 15.19, 15.17, 15.14 and (15.44). □ 



Corollary 15.3. Let 2/3 < 7 < 2. For every e> there is a 6 > such that if x 
constitutes Bianchi IX initial data for l{2.i )- ^2^ ) and 



inf II a: - y II <5 



the 



inf ||$(r, x) 

veA" 

for all T < 0, where $ is the flow of (1^. 



2/11 <e 



Proof. Assuming the contrary, there is an e > and a sequence xi A such that 

inf II $(s,, a;;) -2/11 > e 

yeA 

for some si < 0. Let t; — 0. Since d{Ti, xi) — > and we can assume e is small enough 
that Proposition 14.1 is applicable, there must be an r/ > such that h(si, xi) > rj 



for / large enough, contradicting Theorem 15.1 



□ 



Corollary 15.4. Consider a generic Bianchi IX solution with 2/3 < 7 < 2. Then 
lim {Q + N1N2 + N2N3 + N1N3) ^ 0. 



Proof. If h does not converge to zero, then the conditions of Lemma 15.2 are 
met, since there for a generic solution is an a-limit point on the Kasner circle by 
Proposition 13.1. Corollary 



14.1 



then yields the desired conclusion. □ 



Let A be the set of vacuum type I and II points as in Definition |1.6| . By Corollary 
15.4, a generic type IX solution with 2/3 < 7 < 2 converges to A. 



Corollary 15.5. Let 2/3 < 7 < 2. The closure of J-'ix and the closure ofVix do 
not intersect A. Furthermore, the set of generic Bianchi IX points is open in the 
set of Bianchi IX points. 

Remark. The closure of the Taub type IX points does intersect A. 

Proof. Assume there is a sequence xi G J-ix such that xi x E A. Let t; = 0. 
Observe that then d{xi, ti) — > 0. By Theorem 15.1, there is for each e > and for 
each L a.n I > L such that /i(t, xi) < e for r < t; = 0. By choosing L large enough, 
we can assume n{Ti,Xi) to be arbitrarily small and by choosing e small enough, 
we can assume that Proposition 14.1 is applicable. Consequently, we can assume 



Q,{t,xi) to be as small as we wish for r G (— oo,r;], contradicting the fact that 
r2(T, X;) ^ 1 as T — > —00. The argument for Vix is similar, since the il-coordinate 
of P^{II) is positive. 
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Consider now a generic point x in the set of Bianchi IX points. There is a neigh- 
bourhood of X that does not intersect the Taub points. Let us prove the similar 
statement for Tix. and Viy^. Assume there is a sequence xi £ JFjx such that xi x. 



For each e > there is a T < such that d(T, $(r, x)) < e/2, by Corollary [LK4 . 
By continuity of the flow and the function d, we conclude that for / large enough 
we have d{T,^{T,xi)) < e. Since $(T, a;;) € Tix, we get a contradiction to the 
first part of the lemma. Thus, there is an open neighbourhood of x that does not 
intersect JFjx. The argument for Vix is similar. □ 

Corollary 15.6. Let 2/3 < 7 < 2. The closure of J-yug and the closure ofVyug 
do not intersect A. Furthermore, the generic Bianchi VIIq points are open in the 
set of Bianchi VIIq points. 

Proof. The argument proving the first part is as in the Bianchi IX case, once one has 



checked that analogues of Proposition 14.1 and Theorem 15.1 hold in the Bianchi 



VIIq case. The second part then follows as in the Bianchi IX case, using Proposition 



16. Regularity of the set of non-generic points 



Observe that the constraint (2_^) together with the additional assumption > 
defines a 5-dimensional submanifold of which has a 4-dimensional boundary 
given by the vacuum points. We have the following. 

Theorem 16.1. Let 2/3 < 7 < 2. The sets J-u,J-yUf,, J-ix, 'Pviio o.^T-d 'Pix are 
submanifolds o/M^ of dimensions 1, 2, 3, 1 and 2 respectively. 

We prove this theorem at the end of this section. The idea is as follows. The 
only obstruction to e. g. J-u being a submanifold, is if there is an open set O 
containing F and a sequence Xk G T\i such that x-k — > F , but each x^ has to leave 
O before it can converge to F . If there is such a sequence, we produce a sequence 
yk G such that the distance from to A converges to zero, contradicting 



Lemma 9.1. The argument is similar in the other cases. 

We will need some results from [|o| . The theorem stated below is a special case of 
Theorem 6.2, p. 243. 

Theorem 16.2. In the differential equation 
(16.1) ^'^E^ + G{0 

let G he of class and G{Q) — 0, d^G{Q) — 0. Let E have e > eigenvalues with 
positive real parts, d > eigenvalues with negative real parts and no eigenvalues 



with zero real part. Let — S,it, Co) be the solution of ( 16.1 ) satisfying ^(0, ^o) = Co 
and T* the corresponding map T*{^o) = £,{t,S,o)- Then there exists a map R of a 
neighbourhood of ^ = in -space onto a neighbourhood of the origin in Euclidean 
{u,v) -space, where dim(u) — d and dim(u) = e, such that R is with non- 
vanishing Jacobian and RT^R^^ has the form 

Mt \ ^ / e*^Mo + C/(t, uo,wo) 
vt J \ eJQvo + Vit,UQ,vo) 

U, V and their partial derivatives with respect to uq, vq vanish at (uo,vo) = 0. 
Furthermore V = if vq = and U = if uq = 0. Finally \\e^\\ < 1 and 
lle-Q|| < 1. 



(16.2) 
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Let US begin by considering the local behaviour close to the fixed points. 

Lemma 16.1. Consider the critical point F. There is an open neighbourhood O 
of F in M.^ , and a 1-dimensional submanifold Mu C J^n of O (iTu, such that 



for each x £ O (1 In, either x € Mn, or x will leave O as the flow of (2.1)-(2.o) 
is applied to x in the negative time direction. Similarly, we get a 2-dimensional 
submanifold Afyiio of O RXviiq, and a 3-dimensional submanifold Mix of 
O n Xix with the same properties. Consider the critical point P^{II). We then 
have a similar situation. Give the neighbourhood corresponding to O the name P , 
and use the letter N instead of the letter M to denote the relevant submanifolds. 
Then N^iia has dimension 1 and Nix ho-s dimension 2. 



Proof. Observe that when > 0, we can consider (2.1)-( |2.3D to be an unconstrained 
system of equations in five variables. Using the constraint (2.3) to express fl in 
terms of the other variables, we can ignore fl and consider the first five equations of 
( pj] ) as a set of equations on an open submanifold of , defined by the condition 
f2 > (considering as a function of the other variables). In the Bianchi VIIq 
case, we can consider the system to be unconstrained in four variables. 

Let us first deal with the Bianchi VIIq case. Consider the fixed point P^{II). 
Considering the Bianchi VIIq points with A^i, N2 > Q and A^3 = 0, the linearization 
has one eigenvalue with positive real part and three with negative real part, cf. 



17|. By a suitable translation of the variables, reversal of time, and a suitable 



definition of G and E in ( 16.1 ), we can consi der a solution to ( p.l[ )-(2.£) converging 



to Pi {II) as r ^ —00 as a solution ^ to ( 16.1 ) converging to as t ^ 00. E 



has one eigenvalue with negative real part and three with positive real part, so 



that Theorem 16.2 yields a G^ map i? of a neighbourhood of with non-vanishing 
Jacobian to a neighbourhood of the origin in W^, such that the flow takes the form 
(^6]|) where w G M and v G M^. 



Observe that since f = is a fixed point, there is a neighbourhood of that point such 
that the flow is defined for \t\ < 1. There is also an open bounded ball B centered 
at the origin in (uq, fo)-space such that U and V are defined in a neighbourhood 
N of [—1, 1] X B. Let a = ||e^|| and 1/c = ||e^'5||. For any e > 0, we can choose B 
and then N small enough that the norms of U, V and their partial derivatives with 
respect to u and v are smaller than e in A^. Assume B and N are such for some e 
satisfying 

.rC— 11 — a. 

(16.3) e < min{— —,——}. 



Consider a solution ^ to ( 16.1 ) such that R o ^(t) £ B for all t > T. Let {ut, Vt) = 
R{£,{t)) for t > T. We wish to prove that Vt — 0, and assume therefore that ^ 
for some to > T. We have 

\\vto+n\\ > \\e'^Vto+n-l + V(l,Vt„+n-l,Uto+n-l)\\ > 

1 + C 

> C||ut„+„_i|| - e||-yto+n-l|| > ^:^ll«to+n-l||: 



where we have used ( 16.3 ), the fact that V is zero when vo = 0, and the fact that 



{uf,Vt) remain in B for t>T. Thus, 



||wto+n|| > ( ^-4r- ] ll^*oll> 
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which is irreconcilable with the fact that vt remains bounded. 



If (uto, Vto) & B and vto — 0, ( 16.2 ) yields wto+i = and 



||wto + l|| < (a+ ^-^)ll"tnll = ^Y^ll^tol 

Consequently, all points (m, v) E B with i; = converge to (0, 0) as one applies the 
flow. 

We are now in a position to go backwards in order to obtain the conclusions of 
the lemma. The set R~^{B) will, after suitable operations, including non-unique 
extensions, turn into the set P and R~^{{v = 0} C] B) turns into AViIq- One can 
carry out a similar construction in the Bianchi IX case. Observe that one might 
then get a different P, but by taking the intersection we can assume them to be 
the same. The dimension of A^ix follows from a computation of the eigenvalues. 

The argument concerning the fixed point F is similar. □ 



Proof of Theorem 16.1. Let O, Mu and so on be as in the statement of Lemma 16.1 
Observe that if there is a neighbourhood O C O of F such that Tw n O = Mu n O, 
then .7-11 is a submanifold. The reason is that given any x G Tu, there is a T such 



that <i>(T, x) e O for all t <T. By Lemma [16^, we conclude that $(7", x) e Mu. 
Then there is a neighbourhood O' C O of $(r, x) such that O'nJ^n = O'nMu. We 
thus get, for O' suitably chosen, a map i/j : O' ^ M.^ with inverse, sending 
J-ii n O' to a one dimensional hyperplane. If O' is small enough, we can apply 
<!>(— r, •) to it obtaining a neighbourhood of x. By the invariance of J-u, we have 

$(-T, o') n Tu ^ $(-r, o' n Tn). 

In other words, ^{T,ip{-)) defines coordinates on T, O') straightening out Tu. 
The arguments for the other cases are similar. 

Let us now assume, in order to reach a contradiction, that there is a sequence 
Xk e n O such that Xk ^ F but Xk ^ Mu for all k. If we let O' C O be a 
small enough ball containing F, we can assume that \Ni\' > for z = 1, 2, 3 in O', 



cf. the proof of Lemma 4.2. For k large enough, Xk € O' and applying the flow to 
them we obtain points yk G ^ii H dO' . By choosing a suitable subsequence, we can 
assume that yk converges to a type I point y which is not F. Given e > 0, there 
is a T such that ^{—T, y) is at distance less than e/2 from A. For k large enough, 
r, yk) S Til will then be at distance less than e from A. We get a contradiction 



15.6 and 15.5 



to Lemma 9.1. The arguments for .Fyiif, and J^ix are similar, due to Corollaries 



For T^viio E^iid T'lx, we need to modify the argument. Assume there is a sequence 
Xk G P\nio ^ P such that Xk P^{II), but Xk ^ AViiq for all k. By choosing 
P' C P as a small enough ball, we can assume that |iVj|' > in P' for i = 2,3, 



cf. the proof of Lemma 4.1. For k large enough, Xk S P', and applying the flow to 
them we obtain points yk G Pviio ^ dP'. By choosing a suitable subsequence, we 
can assume that yk converges to a type II point y which is not P^{II). li y ^ Tu^ 
we can apply the same kind of reasoning as before, using Proposition |9.l| to get a 



contradiction to the consequences of Corollary |15.6| . If y G J-n we get, by applying 
the flow to the points yk, a sequence Zk G PviIq converging to F. Applying the 
flow again, as before, we get a contradiction. The Bianchi IX case is similar using 
Corollary ^5^ □ 
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17. Uniform convergence to the attractor 



If X constitutes initial data to (2^)- (2^) at t = 0, then we denote the corresponding 
solution S+(t, x) and so on. 

Proposition 17.1. Let 2/3 < 7 < 2 and let K he a compact set of Bianchi IX 
initial data. Then N1N2N3 converges uniformly to zero on K. That is, for all 
e > there is aT such that 

{NiN2N3)iT,x) < e 

for all T < T and all x £ K . 

Proof. Assume that N1N2N3 does not converge to zero uniformly. Then there is 
an e > 0, a sequence Tk ^ —00 and Xk S K such that 

iNiN2N3)iTk,Xk)>e. 

We may assume, by choosing a convergent subsequence, that X}^ ^ Xi^ as k ^ cxD . 
Because of the monotonicity of {NiN2N3){-,Xk), we conclude that 

iNiN2N3){T,Xk)>e. 

for r G [rfc,0]. Thus 

{NiN2N3){t,x,) = lim {NiN2m){T,Xk) > e 

k — ^00 

for all T < 0. We have a contradiction. □ 

Corollary 17.1. Let 2/3 < 7 < 2 and let K be a compact set of Bianchi IX initial 
data. Then for every e > 0, there is a T such that 



f7 + Ei + < 1 ■ 



for all X £ K and t < T. 



Proof. As before. □ 
Consider 

d = n + N1N2 + N2N3 + N3N1. 

Proposition 17.2. Let K be a compact set of generic Bianchi IX initial data with 
2/3 < 7 < 2. Then d converges uniformly to zero on K. 

Proof. Assume that d does not converge to zero uniformly. Then there is an ry > 0, 
a sequence — s- — 00 and a sequence Xk (z K such that 

(17.1) d{Tk,Xk)>r]. 

We now prove that there is no sequence Sk,^ such that Tk^ < s^^ < and 



Assume there is. By Theorem 15.1, there is no i5 > such that maximum of a;fe„) 



in [Tk,^, Sk„] exceeds d for all n. For S small enough, we can apply Proposition 14.1 



to the interval [t^^, Sk„] to conclude that for some n, Q cannot grow in very much 
in that interval either. We obtain a contradiction to (17.1) for 6 small enough and 
n big enough. 
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Thus there is an e > such that 

d{T,Xk) > e 

for all T G [Tk, 0] and all k. Assume Xk ~* x^. Then 

d(r, a;*) = lim d{T,Xk) > e > 

k — *oo 

for all T < 0. But x^, constitutes generic initial data. □ 

18. Existence of non-special o-limit points on the Kasner circle 

We know that there is an a-limit point on the Kasner circle, but in order to prove 
curvature blow up we wish to prove the existence of a non-special a-limit point on 
the Kasner circle. 

Lemma 18.1. Consider a generic Bianchi IX solution with 2/3 < 7 < 2. If it 

has a special point on the Kasner circle as an a-limit point then it has an infinite 
number of a-limit points on the Kasner circle. 



Proof. By applying the symmetries, we can assume that there is an a-limit point 
on the Kasner circle with = (—1,0). Since the solution is not of Taub 



type, (E+,E_) cannot converge to (—1,0) by Proposition 3.1, Thus there is an 
1 > e > such that for each T there is a r < T such that 1 + E+(t) > e. Let 
Tfc —00 be such that T,+ {Tk) —1. 

Let rj > satisfy rj < e. We wish to prove that there is a non-special a-limit point 
on the Kasner circle with 1 -I- < rj. There is a sequence tk < Tfc such that 
1 -I- I]+(tfe) = rj and E^(ifc) < assuming k is large enough. The condition on 
the derivative is possible to impose due to the fact that 1 -I- E+ eventually has to 
become greater than e. Choosing a suitable subsequence of {^fcj, we get an a-limit 



point which has to be a vacuum type I or II point by Corollary 15.4. If it is of type 
I, we get an a-limit point on the Kasner circle with 1 + = rj and we are done. 
The a-limit point cannot have iVi > 0, because of the condition on the derivative, 



cf. the proof of Proposition 5.1. If it is of type II with N2 or greater than 
zero, we can apply the flow to get a type II solution, call it x, of a-limit points 
to the original solution. Since a type II solution with A^2 or greater than zero 
satisfies E^ < 0, the w-limit point y oi x must have 1 -|- E_|_ < rj. By Proposition 



5.1 y G /C2 U /C3, so that it is non-special. 

Let < 7/1 < e. As above, we can then construct a non-special a-limit point xi 
on the Kasner circle with coordinate S+^i such that 1 -I- < 771. Assume 
we have constructed non-special a-limit points Xi on the Kasner circle, i = 1, ...,777 
with S+ coordinates S+^i satisfying S+_i < Let < rjm+i < 1 -|- S+_,„. 

Then by the above we can construct a non-special a-limit point Xm+i on the Kasner 
circle with E_|_ coordinate E+.m+i, satisfying E+^m+i < E+,m- Thus the solution 
has an infinite number of a-limit points on the Kasner circle. □ 

Corollary 18.1. A generic Bianchi IX solution with 2/3 < 7 < 2 has at least three 
non-special a-limit points on the Kasner circle. Furthermore, no Ni converges to 
zero. 



Proof. Assume first that the solution has a special a-limit point on the Kasner 
circle. By Lemma 18.1, the first part of the lemma follows. By the proof of 
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Lemma 18.1, there is a non-special a-limit point on the Kasner circle with 
coordinate arbitrarily close to —1, say that it belongs to /C2- Repeated application 
of Proposition 6.1 then gives a-limit points first in /C3, and after enough iterates, 
either an a-limit point in /Ci, or a special a-limit point on the Kasner circle with 
1/2. If the latter case occurs, a similar argument to the proof of Lemma 
18.1 yields an a-limit point on JCi. By Proposition |6.l|, we conclude that there are 



a-limit points with A^i > 0, with A^2 > and with A^3 > 0. 

Assume that there is no special a-limit point on the Kasner circle. Repeated 
application of the Kasner map yields a-limit points in ICi, i = 1,2,3, and the 
conclusions of the lemma follow as in the previous situation. □ 



19. Conclusions 

Let us first state the conclusions concerning the asymptotics of solutions to the 
equations of Wainwright and Hsu. We begin with the stiff fluid case. 



Theorem 19.1. Consider a solution to (2.1)-(2.o) with 7 = 2 and Q > 0. Then 
the solution converges to a type I point with -|- < 1. For the Bianchi types 
other than I, we have the following additional restrictions. 

1. // the solution is of type II with Ni > 0, then S+ < 1/2. 

2. For a type VIq or VIIq with N2 and non-zero, then ± -v/SS- > — 1. 

3. // the solution is of type VIII or IX, then ± \/35]_ > —1 and < 1/2. 

Remark. Figure [s] illustrates the restriction on the shear variables. The types 
depicted are 1, 11, VIo and VIIq, and Vlll and IX, counting from top left to bottom 
right. 



Proof. The theorem follows from Propositions 7.1 and 7.2. □ 



Consider now the case 2/3 < 7 < 2. Let A be the closure of the type II vacuum 
points. 

Theorem 19.2. Consider a generic Bianchi IX solution x with 2/3 < 7 < 2. Then 
it converges to the closure of the set of vacuum type II points, that is 

lim inf ||x(t) — y|| = 

r— > — 00 y^A 

where \\ ■ \\ is the Euclidean norm on M^. Furthermore, there are at least three 
non-special a-limit points on the Kasner circle. 



Remark. One can start out arbitrarily close to this set without converging to it, cf. 



Proposition 11.1 



Proof. The first part follows from Corollary 15.4 and the second part follows from 
Corollary pOl □ 



Proof of Theorems 1.1 and l.i.. Let {M,g) be the Lorentz manifold obtained in 



Lemma 21.2 with topology I x G. It is globally hyperbolic by Lemma 21.4 



If the initial data satisfy tr^fc — for a development not of type IX, then it is 
causally geodesically complete and satisfies ^ = for the entire development, by 



Lemma 21.5 and Lemma 21.8. The first part of Theorem 1.1 follows 
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Figure 8. The points to which the shear variables may converge 



for a stiff fluid. 



Cons ider initial data o f type I, II, VIo, VIIo or VIII such that tr^fc 0. By Lemma 
21.5 and Lemma 21. 8| , we may then time orient the development so that it is future 
causally geodesically co mple te and past causally geodesically incomplet e, and the 
second part of Theorem 1.1 follows. The third part follows from Lemma 21.8. 



Consider an inextendible future directed causal geodesic in the above dev elopm ent. 
Since each hypersurface {v} x G is a Cauchy hypersurface by Lemma 21.4 , the 
causal curve exhausts the interval /. 

1. If the solution is not of type IX, then the solution to ( 21.4 )-( ^L9 ), which is used 
in constructing the class A development, corresponds to a solution to (2.1)-(2.3), 
because of Lemma 21.5. Furthermore, i — > t_ corresponds to r ^ — oo, because of 



Lemma 22.4 



a. In all the stiff fluid cases, the solution to (2.1)-(2.3) converges to a non-vacuum 
type I point by Theorem 19.1, so that Lemma 22.1 and Lemma 22.3 yield the 



desired conclusions in that case. 

b. Type I, II and VIIq with 1 < 7 < 2. That the Kretschmann scalar is unbounded 



Proposition 10.2, Lemma 22.1 and and Lemma 22.2 



in the cases stated in Theorem 1.2 follows from Proposition S.l, Proposition 9.1 



c. Non-vacuum solutions which are not of type IX. Then RapR"^ is unbounded 
xising Lemma 22.3. 



THE BIANCHI IX ATTRACTOR 



67 



2. If the solution is of type IX, then half of a solution to (21.4)-( 21.9| ) corresponds 
to a Bianchi IX solution to (2.1)-(2.3); because of Lemma 21. £. By Lemma ^2.5 , 
t ^ t± corresponds to r 



-oo. 



a. In the stiff fluid case, we get the desired statement as before. 

b. If 1 < 7 < 2, we get the d esired conclusions , con cerning blow up of the 
Kretschmann scalar, from Corollary 18.1 , Proposition |ll.l| , Lemma 22.1 and Lemma 



22.2 



c. Non-vacuum solutions. Then RajsR"^ is unbounded using Lemma 22.3 



Let us now prove that the development is inextendible in the relevant cases. Assume 
there is a connected Lorentz manifold (M, g) of the same dimension, and a map 
i : M M which is an isometry onto its image, with i(A'I) ^ M . Then there is 
a p S M — i(Af ) and a timelike geodesic 7 : [a, 6] ^ M such that 7([a, 6)) C i(M) 
and 7(6) = p. Since 7|[a,6) can be considered to be a future or past inextendible 
timelike geodesic in M, either it has infinite length or a curvature invariant blows 
up along it, by the above arguments. Both possibilities lead to a contradiction. 
Theorem 1.2 follows. □ 



20. Asymptotically velocity term dominated behaviour near the 

singularity 



In this section, we consider the asymptotic behaviour of Bianchi VIII and IX stiff 
fluid solutions from another point of view. We wish to compare our results with 
|2j, a paper which deals with analytic solutions of Einstein's equations coupled to a 
scalar field or a stiff fluid. In Andersson and Rendall prove that given a certain 
kind of solution to the so called velocity dominated system, there is a unique solution 
of Einstein's equations coupled to a stiff fluid approaching the velocity dominated 
solution asymptotically. We will be more specific concerning the details below. The 
question which arises is to what extent it is natural to assume that a solution has 
the asymptotic behaviour they prescribe. We show here that all Bianchi VIII and 
IX stiff fluid solutions exhibit such asymptotic behaviour. 

In order to speak about velocity term dominance, we need to have a foliation. In our 
case, there is a natural foliation given by the spatial hypersurfaces of homogeneity. 
Relative to this foliation, we can express the metric as in ( ^1.14 ) according to 
Lemma 21.2. In what follows, we will use the frame e\ appearing in Lemma ^1.2 , 
and Latin indices will refer to this frame. Let g be the Riemannian metric, and fc 
the second fundamental form of the spatial hypersurfaces of homogeneity, so that 



(20.1) 



9ij ^gie'i: 



where g is as in (21.14). The constraint equations in our situation are 
(20.2) R - kijk'^ + [ixkf = 2fj. 



(20.3) 



V,(trfc) 



0, 



which are the same as ( 21.8 ) and ( 21. 5| ) respectively. The evolution equations are 



(20.4) 
(20.5) 



dtk\ 



= R\ 4 



(trfc)fc\. 
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The evolution equation for the matter is 



(20.6) 



dt^ — 2(trfc)/i. 



We wish to compare solutions to these equations with solutions to the so called 
velocity dominated system. This system also consists of constraints and evolution 
equations, and we will denote the velocity dominated solution with a left superscript 
zero. The constraints are 



(20.7) 
(20.8) 



- o%Ofc'^' + (tr0fc)2 = 2V 
^V*(°%) -°Vj(tr"fc) = 0. 



The evolution equations are 

(20.9) = -2°% 

(20.10) at°fc^ = (tr°fc)°fc^, 
and the matter equation is 

(20.11) = 2(tr°fc)V- 

We raise and lower indices of the velocity dominated system with the velocity 
dominated metric. In [ 2|, Andersson and Rendall prove that given an analytic 
solution to (20.7)-(20.11) on S* x (0, oo) such that ttr^fc = —1, and such that the 
eigenvalues of —t^k'^j are positive, there is a unique analytic solution to ( ^0.2 )-(20.6) 
asymptotic, in a suitable sense, to the solution of the velocity dominated system. In 
fact, they prove this statement in a more general setting than the one given above. 
We have specialized to our situation. Observe the condition on the eigenvalues of 
—t^k^j. Our goal is to prove that this is a natural condition in the Bianchi VIII 
and IX cases. 

Theorem 20.1. Consider a Bianchi VIII or IX stiff fluid development as in Lemma 



21. H: with fiQ > 0. Choose time coordinate so that i_ = 0. Then there is a solution 
to { 20.1 )-{ 20.11 } such that ttr°fc — —\, the eigenvalues of —t'^k^ ^ are positive, and 



the following estimates hold 



1. Vsij- =<5', +o(i"'0 

2. fc^j = "fc^^. +o(i-i+"'0 

3. /i = V + o(i-2+/3i), 



where ^ and /3i are positive real numbers. 

Remark. In |^ two more estimates occur. They are not included here as they are 
replaced by equalities in our situation. Observe that the difhculties encountered in 
j2| concerning the non-diagonal terms of fc*^ disappear in the present situation. 



Proof. Below we will use the results of Lemma 21.2 and its proof implicitly. When 



we speak of aij, 9, riij and n, we will refer to the solution of (21.4)-(21.9) and 



the indices of these objects should not be understood in terms of evaluation on a 
frame. Since 9ij and so on are all diagonal, we will sometimes write 9i etc instead, 
denoting diagonal component i. The re a re two relevant frames: e[ and — a^e^. 
The latter frame yields riij through (1.7). When we speak of fc*^-, Rij and so on, 
we will always refer to the frame e-. We have 
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(no summation on i). The metric is given by ( ^0.1 ) above. Let us choose 



(20.12) 

let °9 = 06*1 



and 



^3 



-2 r 



(no summation on i). Because of (20.12), equation ( 20.^ ) will be satisfied since it is 
a statement concerning the commutation of "fc*^- and riij. The existence interval for 
the solution to Einstein's equations is (0, t+) by our conventions, and since we wish 
to have ttr^fc = —1 we need to define '^9{t) = l/t. Observe that is constant 

in time, and that 9i/9 converges to a positive value as < — *■ 0; this is a consequence 
of Theorem 19.1 and the definition ( ^1.11 ) of the variables S+ and E_. Choose °6'i 
so that ^9i/'^9 coincides with the limit oi9i/9. Similarly is constant, ^,/9^ 

converges to a positive value, and we choose ^^/^9^ to be the limit. Since R/9^ 
is a polynomial in the Ni and the Ni converge to zero by Theorem 19.1, equation 
( pO^ ) will be fulfilled. By our choices, ( p0T0| ) and ( pOTl| ) will also be fulfilled. We 
will specify the initial value of ^a^ later on, and then define "a^ by demanding that 



( pO.gj ) holds. 

It will be of interest to estimate terms of the form i?' ^ / 9^ . These terms are quadratic 
polynomials in the Ni. By abuse of notation, we will write Nilr) when we wish 
to evaluate Ni in the Wainwright-Hsu time ( 21.10| ) a nd N i{t) when we wish to 
evaluate in the time used in this theorem. By Theorem 19.1 , there is an e > and 
a To such that 

\N,{t)\ <exp(eT) 

for all T < To. We wish to rewrite this estimate in terms of t. Let us begin with 
( 21.12| ). Since we can assume that g < 3 for t < To we get 



r) < exp[-4(T-To)]%o), 



so that for ti, t < To we get, using (21.10) 



t{T) - tin) = ^ ^ds> ^^(exp[4(T - To)] - exp[4(Ti - To)]). 



Letting ti go to —oo and observing that t(~oo) — 0, cf. Lemma 22.4 and Lemma 
22.5, we get for some constant c 



so that 



e^" < ct(T), 



N,{t) < exp(eT(t)) < Ct"^ 



for some positive number rj. Consequently expressions such as R''j/9^ and R/9^ 
satisfy similar bounds. 

Let us now prove the estimates formulated in the statement of the theorem. Observe 
that for t small enough, we have 



^9 = trfc(t) = -( 



R 



l]ds)-\ 
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since the singularity is at i = and trfc must become unbounded at the singularity, 
cf. Lemma ^2-4 ^2.5| and ( |21.12D . Thus we get 

(20.13) e-^e = -f^ + ^1^'^"' ^ or'+"^) 

for some ?7i > 0. In order to make the estimates concerning fc*^-, we need only 
consider 6i and ^Oi. We have 



dti- 



iR — W 



Oa. 0/ 



with no summation on the i in This computation, together with the estimates 
above and the fact that Oi/9 — °0i/°6' converges to zero, yields the estimate 

(20.14) 

for some 772 > 0. However, 
(20.15) + 

Combining ( 20.13| ), (20.14) and ( 20.15| ), we get estimate 2 of the theorem. Similarly, 
we have 

*^5)2 002' *02 dis- 
integrating, using the fact that converges to "/.t/^^^, we get 



(20.16) 

where 773 > 0. Using 



Q2 002 



o(f'^) 



_M 

02 002 



V ^i - V 



02 002 02 ' 

( 20.13| ) and ( 20.16D , we get estimate 3 of the theorem. Finally, we need to specify 
the initial value of and prove estimate 1. Since 

(no summation on i) and similarly for '^Ui, we get 



„ "I _ U-i ,g 

Ct^— - 75 — [ 



By our estimates on ^9i — 9i, we see that this implies that ai/°ai converges ast ^ 0. 
Choose the value of "a^ at one point in time so that this limit is 1. We thus get, 
using estimate 2 of the theorem. 

Estimate 1 of the theorem now follows from this estimate and the fact that 



^1 S' 



The theorem follows. □ 
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21. Appendix 

The goal of this appendix is to relate the asymptotic behaviour of solutions to the 



ODE (2T)-(2^) to the behaviour of the spacetime in the incomplete directions of 



inextendible causal curves. We proceed as follows. 

1. First, we formulate Einstein's equations as an ODE, assuming that the space- 



time has a given structure ( 21.1 ). The first formulation is due to Ellis and 
MacCallum. We also relate this formulation to the one by Wainwight and 
Hsu. 



2. Given initial data as in Definition 1.1, we then show how to construct a 



Lorentz manifold as in (21.1), satisfying Einstein's equations and with initial 
data as specified, using the equations of Ellis and MacCallum. We also prove 
some properties of this development such as Global hyperbolicity and answer 
some questions concerning causal geodesic completeness. 



3. Finally, we relate the asymptotic behaviour of solutions to (2.1)-(2.3) to the 
question of curvature blow up in the development obtained by the above 
procedure. 

We consider a special class of spatially homogeneous four dimensional spacetimes 
of the form 

(21.1) (M,5) - (/ X G, -dt^ + x.^m ® e), 

where / is an open interval, G is a Lie group of class A, Xij is a smooth positive 
definite matrix and the are the duals of a left invariant basis on G. The stress 
energy tensor is assumed to be given by 

(21.2) T = ^idt^ +p{g + dt^), 

where p = (7 — l)/i. Below, Latin indices will be raised and lowered by Sij. 



Consider a four dimensional {M,g) as in (|21.l| ) with G of class A. In order to 
define the different variables, we specify a suitable orthonormal basis. Let eg = dt 
and Ci = a^-' Zj, i=l,2,3, be an orthonormal basis, where a is a C°° matrix valued 
function of t and the Zi are the duals of . 

By the following argument, we can assume that < Veo&ii^j 0. Let the matrix 
valued function A satisfy eo(A) + AB = 0, A(0) — Id where Bij =< Vf^^Ci^ej > 
and Id is the 3x3 identity matrix. Then A is smooth and 50(3) valued and if 
e- = A/ej, then < Veoe-,e^ >= 0. 

Let 

(21.3) e{X,Y)=<Vxeo,Y>, 

ep) and [e/3,e^] = Ip^ea where Greek indices run from to 3. The 
objects 0af3 and 7^^ will be viewed as smooth functions from / to sonic suitable 
M*^, and our variables will be defined in terms of them. 

Observe that [Zi, cq] ~ 0. The span the tangent space of G, and < [eo, e^], eo 

0. We get 6qo = Oot = Q and Oa/i symmetric. We also have 7° — — and 

7oj ~ We let n be defined as in ( |l.7| ) and 
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where we by abuse of notation have written tr(0) as 9. 

We express Einstein's equations in terms of n, a and 9. The Jacobi identities for 
Ca yield 



—9nji — 0. 
3 ^ 



(21.4) eo(ny) 

The Oi-components of the Einstein equations are equivalent to 

(21.5) cr^nfcj - n^'^akj = 0. 

Letting bij — 2n^'^nkj — tr{n)nij and Sij — bij — itr(6)(5,j , the trace free part of the 
ij equations are 

(21.6) eo(crij) + 9aij + Sij = 0. 
The 00-component yields the Raychaudhuri equation 

1 



(21.7) 



eo(^) 



e,j9'^ + -{3-1-2)^1^0, 



and using this together with the trace of the ij-equations yields a constraint 
(21.8) (Tija'^ + {nijn'^ - itr(n)2) + 2fi = ^9^ . 



Equations ( 21.4 )-( 21.8 ) are special cases of equations given in Ellis and MacCallum 
|8|. A t a point to, we may diagonalize n and a simultaneously since they commute 
(21.5). Rotating Bq by the corresponding element of 5*0(3) yields upon going 
through the definitions that the new n and a are diagonal at to. Collect the off- 
diagonal terms of n and a in one vector v. By (21.4) and (21.6), there is a time 
dependent matrix C such that v = Cv so that v{t) — for all t, since v{to) = 0. 
Since the rotation was time independent, < \I ^^Ci, ej >~ holds in the new basis. 

The fact that T is divergence free yields 

(21.9) eo{tJ.)+-/9fi = 0. 



Introduce, as in Wainwright and Hsu |17|, 

^(Jij/9 
Nij = nij/9 

n = 3fi/9^ 

and define a new time coordinate r, independent of time orientation, satisfying 

eft _ 3 
d^~9' 

For Bianchi IX developments, we only consider the part of spacetime where 9 is 
strictly positive or strictly negative. Let 



(21.10) 



(21.11) 



-33, 



and E_ = ^{^22 - S33). 



If we let Ni be the diagonal elem ents of Nij, equations ( 21.4 ) and ( 21.6 ) turn into 
( ^^ ) with definit ions as in (2_^), except for t he expressio n fo r il' . It can however 
be derive d from ( 21.9 ). The constraint ( 21.8 ) turns into (2_^). The Raychaudhuri 
equation ( 21.7 ) takes the form 

(21.12) 9' = -{l + q)9. 
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Before using the equations of Ellis and MacCallum to construct a development, it 
is convenient to know that one can make some simplifying assumptions concerning 
the choice of basis. The next lemma fulfills this objective, and also proves the 
classification of the class A Lie algebras mentioned in the introduction. 

Lemma 21.1. Table^ constitutes a classification of the class A Lie algebras. Con- 
sider an arbitrary basis {ci} of the Lie algebra. Then by applying an orthogonal 
matrix to it, we can construct a basis {e^} such that the corresponding n' defined 



by (LI) is diagonal, with diagonal elements of one of the types given in Table |]. 



Proof. Let be a basis for the Lie algebra and n be defined as in (1.7). If we 



{A- 



'cj, then n transforms to 



change the basis according to e- 

(21.13) 7i' = (det^)" 

Since n is symmetric, we assume from here on that the basis is such that it is 
diagonal. The matrix A = diag(l 1—1) changes the sign of n. A suitable 
orthogonal matrix performs even permutations of the diagonal. T he num ber of non- 
zero elements on the diagonal is invariant under transformations (^1.13[) taking one 
diagonal matrix to another. If A = (oij) and the diagonal matrix n' is constructed 
as in (21.13), we have n'j.^. — {det A)~^Y!,^=i'^ik''^iii so that if all the diagonal 
elements of n have the same sign, the same is true for n' . The statements of the 
lemma follow. □ 



We now prove that if we begin with initial data as in Definition 1.1, we get a 
development as in Definition 1.4 of the form ( ^1.1 ), with certain properties. 



Lemma 21.2. Fix 2/3 < 7 < 2. Let G, g, k and fiQ be initial data as in Definition 
Then there is an orthonormal basis e' i = 1,2,3 of the Lie algebra such that 



defined by (Uj) and kij ~ k{e[, e'j) are diagonal and n'^^ is of one of the forms 



given in Table IJ. Let 



1 



61(0) = -tvgk, cry(O) = -k{e'i,e'j) + -6'(0)%, 71^(0) = n'^^ and /i(0) = a^o- 

Solve (21.4), ( 21. (\ ), (21."}) and ( 21.1^ with these conditions as initial data to obtain 
n, a, and fi, and let I be the corresponding existence interval. Then there are 
smooth functions ai : I (0, 00) i — 1,2, 3, with ai{0) = 1, such that 



(21.14) 



9 



3 



dt' + }_^ar^it)C®C, 



where is the dual of e[, satisfies Einstein's equations (1.5) on M = I x G, with 
T as in (1.1) with u — cq, /i as above and p — (7 — l)/i. Furthermore, 



< Ve,eo, e, > = 



-eSi 

3 ■ 



where V is the Levi-Civita connection of g and Ci = aie[, if we consider the left hand 
side to be a function oft. Consequently, the induced metric and second fundamental 
form on {0} x G are g and k, and we have a development satisfying the conditions 
of Definition l.J). 



Proof. Let e^, i = 1, 2, 3 be a left invariant orthonormal basis. We can assume 
the corresponding n' to be of one of the forms given in Table |l| by Lemma [21.1 



The content of (1.5) is that kij — fc(e^,e') and n' are to commute. We may 
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thus also assume fc^ to be diagonal without changing the earlier conditions of the 
construction. If we let n(0) — n', 0(0) = — tr^/c, (7^(0) — —kij + 6Sij/3 and 
l40)= ^J■o, t hen ( pT^ ) is the same as ( |21.8| ). Let n, cr, 6* and /i s atisfy ( |21.6D 



( ^1.7 ) and ( 21. 9| ) with initial values as specified above. Since ( 21. 8| ) is satisfied at 



0, it is satisfied for all times. For reasons given in connection with (21.8), n and 



cr will remain diagonal so that (21.5) will always hold. Let rii and at denote the 
diagonal elements of n and a respectively. 

How are we to define the in the statement of the lemma? The n obtained from 



Ci by (1.7) should coincide with n. This leads us to the following definitions. Let 
/i(0) = 1 and fi/ fi = 2ai — 9/3. Let — (nj^i/j)^/^ and define = 0^6-. Then 
h associated to Ci equals n. We complete the basis by letting eo — 9*. Define 
a metric < •, • > on M by demanding Ca to be orthonormal with eo timelike 
and Ci spacelike, and let V be the associated Levi-Civita connection. Compute 
< Veoei,ej >= 0. Ee{X,Y) \7xeo,Y > and 0^,^ = 0{e^,e^), then §00 = = 
9oi — 0. Furthermore, 

aj 

(no summation over j) so that 9ij is diagonal and tr6 = 6. Finally, 

— CTii = —da + 3^ ~ 

The lemma follows by considering the derivation of the equations of Ellis and Mac- 
Callum. □ 



Definition 21.1. A development as in Lemma 21.2 will be called a class A devel- 
opment. We will also assign a type to such a development according to the type of 
the initial data. 

The next thing to prove is that each = {v} x G is a Cauchy surface, but first 
we need a lemma. 

Lemma 21.3. Let p be a left invariant Riemannian metric on a Lie group G. Then 
p is geodesically complete. 

Proof. Assume 7 : (t_, i+) — > G is a geodesic satisfying p{-^' , 7') = 1, with i+ < 00. 
There is a 5 > such that every geodesic A satisfying A(0) = e, the identity element 
of G, and A'(0) = v with p(u, w) < 1 is defined on (— 5, 5). If Lh : G ^ G is defined 
by Lh{hi) = hhi, then L^ is by definition an isometry. Let to G (t_,t+) satisfy 
t+^to < (5/2. Let V e TeG be the vector corresponding to 7'(to) under the isometry 
Let A be a geodesic with A(0) = e and A'(0) = v. Then L^(^t^-^)oX is a geodesic 
extending 7. □ 

Let us be precise concerning the concept Cauchy surface. 

Definition 21.2. Consider a time oriented Lorentz manifold {M,g). Let / be an 
interval in R and 7 : / — > M be a continuous map which is smooth except for a finite 
number of points. We say that 7 is a future directed causal, timelike or null curve if 
at each t E L where 7 is differentiable, j'it) is a future oriented causal, timelike or 
null vector respectively. We define past directed curves similarly. A causal curve is 
a curve which is either a future directed causal curve or a past directed causal curve 
and similarly for timelike and null curves. If there is a curve A : /i — ^ Af such that 
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7(1) is properly contained in A(/i), then 7 is said to be extendible^ otherwise it is 
called inextendible. A subset 5* C M is called a Cauchy surface if it is intersected 
exactly once by every inextendible causal curve. A Lorentz manifold as above which 
admits a Cauchy surface is said to be Globally hyperbolic. 



Lemma 21.4. For a class A development, each My — {v} xG is a Gauchy surface. 



Proof. The metric is given by ( ^1.14 ). A causal curve cannot intersect twice 



since the t-component of such a curve must be strictly monotone. Assume that 
7 : (s_,s+) — > M is an inextendible causal curve that never intersects M.^. Let 
t : M ^ / be defined by t[{s, h)] — s. Let sq S (s-, s+) and assume that ^(7(30)) = 
ti < V and that < 7', 9* >< where it is defined. Thus t{j{s)) increases with s 
and t(7([so, s+))) C Since we have uniform bounds on from below and 

above on [ti, v] and the curve is causal, we get 

3 

(21.15) {Y.chry^'<-c<y,eo> 

i=l 

on that interval, with C > 0. Since 

(21.16) / -<7',eo>ds=/ —r-^ds<v-ti, 

the curve 7|[so,s+)i projected to G, will have finite length in the metric p on G 
defined by making an orthonormal basis. Since p is a left invariant metric on 



a Lie group, it is complete by Lemma 21.3, and sets closed and bounded in the 



corresponding topological metric must be compact. Adding the above observations, 
we conclude that 7([so,s+)) is contained in a compact set, and thus there is a 
sequence Sk S [so,s+) with Sk s+ such that 7(5^) converges. Since ^(7(3)) is 



monotone and bounded it converges. Using (21.15) and an analogue of (21.16), 
we conclude that 7 has to converge as s ^ s+. Consequently, 7 is extendible 
contradicting our assumption. By this and similar arguments covering the other 
cases, we conclude that My is a Cauchy surface for each v G (i_,t+). □ 

Before we turn to the questions concerning causal geodesic completeness, let us 
consider the evolution of 9 for solutions to the equations of Ellis and MacCallum. 
This is relevant also for the definition of the variables of Wainwright and Hsu, since 
there one divides by 9. We first consider developments as in Lemma 21. 2| which are 
not of type IX. 

Lemma 21.5. Gonsider class A developments which are not of type IX. Let the 
existence interval be I = (t_,i+). Then there are two possibilities. 

1. 9 ^ for the entire development. We then time orient the manifold so that 
9 > Q. With this time orientation, t+ = 00. 

2. 9 — 0, aij — and fi — for the entire development. Furthermore, ntj is 
constant and diagonal and two of the diagonal components are equal and the 
third is zero. The only Bianchi types which admit this possibility are thus type 
I and type VIIq. Furthermore I = (— cx),oo). 
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Proof. Sin ce n^, is diagonal, see the proof of Lemma 21.2, we can formulate the 
constraint ( 21.8| ) as 



1 



ana ■ 



+ -[ni + (ns - n^y - 2ni(n2 + ng)] + 2/i 



where the rii wee the diagonal components of n^ . Considering Table ^ we see that, 
excepting type IX, the expression in the rii is always non-negative. Thus we deduce 
the inequality 

2 



(21.17) 



a^ja'^ +2fi< -e^ 



Combining it with ( |21.7| ), we get |eo(6')| < 9"^, using the fact that 2/3 < 7 < 2. 
Consequently, if 6 is zero once, it is always zero. Time orient the developments 
with 6* 7^ so that 9 > 0. 



Consider the possibility 9 — 0. Equa tion ( |21.7| ) then implies cr.y — and ^ = 



since 7 > 2/3. Equations ( 21.8 ) and ( 21.6 ) then imply bij — 0, and ( 21.4 ) implies 



constant. All the statements except the the fact that t-i- = 00 in the 9 > case 
follow from the above. 

Observe that 9 decre ases in magnitude with time, so that it is bounded to the 
future. By the ( 21.17 ), the same is true of aij and /i. Using ( 21.4 ), we get control 
of riij and conclude that the solution may not blow up in finite time. We must thus 
have = 00. □ 

By a theorem of Lin and Wald |0, Bianchi IX developments recoUapse. 

Lemma 21.6. Consider a Bianchi IX class A development with 1 < 7 < 2 and 
I — (i_ , i+). Then there is a to d I such that 9 > in (i_, to) and 9 < in {to, t^). 

Proof. Let us begin by proving that 9 can be zero at most once. If 9(ti) = 0, 
i = 1,2 and ti < t2, then 6* = in (ti,t2) since it is monotone by (21.7). Thus 
(21.7) implies (Ty = = /i in (ti,t2) as well. Combining this fact with (21.8) and 
(21.(:), we get bij = 0, which is impossible for a Bianchi IX solution. Assume 9 is 
never zero. By a suitable choice of time orientation, we can assume that 9 > on 
I. Let us prove that t+ = 00. Since 9 is decreasing on Ii = [0, t+) and non-negative 
on / it is bounded on Ii. By (21.4), ?ii7i2'T-3 decreases so that it is bounded on Ii. 



By an argument similar to the proof of Lemma 3.3, one can combine this bound 
with ( 21. 8| ) to conclude that aij and fi are bounded on Ii. By (21.4), we conclude 



that Hij cannot grow faster than exponentially. Consequently, the future existence 
interval must be infinite, that is i+ = 00, since / was the maximal existence interval 
and solutions cannot blow up in finite time. In order to use the arguments of Lin 
and Wald, we define 



a,{s)ds + f]'^ , a{t) 



1 



9{s)ds + ao, 



where 2/3° - ao = ln(n,(0)) and f^? = 0- Then 

Hi = exp(2/3i - a). 

Let p — ^/Sir and Pi — p/Sir = (7 — l)^/87r, i ~ 1,2,3. Equations ( 21. 8| ) and 
(21.7) then imply equations (1.4) and (1.5) of |l^, and equations (1.6) and (1.7) 
of [IJI follow from ( ^1.6 ). We have thus constructed a solution to (1.4)-(1.7) of 
14| on an interval [0,oo) with da/dt > 0. Lin and Wald prove in their paper [|3 
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that this assumption leads to a contradiction, if one assumes that \Pi\ < p and 
Pi+ P2+ > 0. However, these conditions are fulfilled in our situation, assuming 
1 < 7 < 2. In other words, there is a zero and since is decreasing it must be 
positive before the zero and negative after it. The lemma follows. □ 

The lemma concerning causal geodesic completeness will build on the following 
estimate. 



Lemma 21.7. Consider a class A development. Let 7 : (s_, s+) — > 
directed inextendible causal geodesic, and 

(21.18) /,(,s)-<7'(s),e,|^(,) >. 

If ~ for the entire development, then /o is constant. Otherwise, 



M he a future 



(21.19) 



V2, 



Remark. We consider functions of t as functions of s by evaluating them at i{j{s)), 



where t is the function defined in Lemma 21.4 



Proof. Compute, using the proof of Lemma 21. 2| , 

3 

=<7'(s),Vy(,)eo >^Y,dkfi, 



dfo 
ds 



k=l 



where 9k are the diagonal elements of 9ij. If 9 = for the entire devel opment, then 
9k = for the entire development by Lemma 21.5 and Lemma 21.6, so that /o is 
constant. Compute, using Raychaudhuri's equation (21.7), 

33 3 



fc=l k=l fe=l 



i(37-2)/x/o^ 



where ak are the diagonal elements of . Estimate 

1/2/3 \ 3 

E-0 E/- 

\k=l / fe=l 

using the tracelessness of aij. By making a division into the three cases J2k=i'^k — 
6*2/3, 6*73 ^ELi CTfe < 26*73 and 26lV3 < XlLi '^h and using the causality of 7 
we deduce ( |21.19| ). □ 

Lemma 21.8. Consider a class A development with existence interval I = (t_,t+). 
There are three possibilities. 

1. 9 = for the entire development, in which case the development is causally 
geodesically complete. 

2. The development is not of type IX and 9 > 0. Then all inextendible causal 
geodesies are future complete and past incomplete. Furthermore, t_ > ~oo 
and t^ ~ 00. 

3. // the development is of type IX with 1 < 7 < 2, then all inextendible causal 
geodesies are past and future incomplete. We also have i_ > ~oo andt^ < 00. 

Proof. Let 7 : {s^, s+) M be a future directed inextendible causal geodesic and 
fi, be defined as in (|21.18). Let furthermore / = (t_,t+) be the existence interval 



mentioned in Lemma 21.2. Since every My, u e / is a Cauchy surface by Lemma 
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|21.4| , t(7(s)) must cover the interval I as s runs through {s-,S-^.). Furthermore, 
t(7(s)) is monotone increasing so that 

(21.20) i{j{s)) ^ t± as s ^ s±. 

Let So G (s-jS-i-) and compute 



(21.21) 



'fo{u)du = t(7(s)) - i(7(so)). 



Consider the case 6 — for the entir e dev elopment. B y Lem ma 21.7, fp is then con- 
stant, and / — (— oo, cx)) by Lemma 21.5. Equations ( |21.21 ) and (21.20) then prove 
that we must have (s_,s+) = (— cx),cx)). Thus, all inextendible causal geodesies 
must be complete. 

Assume that the development is not of type IX and that 9 > 0. S ince f pO is negative 
on [so,s+), its absolute value is bounded on that interval by (|21.19| ). If 5+ were 
would be bounded from below by a positive constant on [so, s+), since 

if I < < ce 



finite, 



on that interval for some C > 0, cf. (21.17) and the observations following that 
equation. Since fpO is bounded, we then deduce that /□ is bounded on [sq,s+). 
But then ( |21.20| ) and ( ^1.21| ) cannot both hold, since t+ = oo by Lemma ^1.5| 



Thus, s+ = oo and all inexte ndible causal geodesies must future complete. Since 



fo9 is negative on (s_,s+), ( 21.19 ) proves that this expression must blow up in 
finite s-time going backward, so that s_ > —oo. Since the curve 7(5) = (s,e) is an 
inextendible timelike geodesic, we conclude that t_ > —00. 

Consider the Bianchi IX case. By Lemma 21.4 and pl.t: , we conclude the existence 
of an So G (s_,s+) such that fo9 is negative on (s_,So) and positive on (so,s_|_). 
By ( 21.19D , fo9 must blow up a finite s-time before so, and a finite s-time after Sq. 
Every inextendible causal geodesic is thus future and past incomplete. We conclude 
t_ > —00 and < 00. □ 



22. Appendix 



In this appendix, we consider the curvature expressions. According to p9|] , p. 40, 
the Weyl tensor Cap-yS is defined by 

RaPtS = Caf3-fS + (5q[7-R<5]/3 " gf3[-yRs]a) " -^R9a[-y9S]l3 , 

where the bar in gap and so on indicates that we are dealing with spacetime objects 
as opposed to objects on a spatial hypersurface. Using this relation and the fact 
that our spacetime satisfies (1.3), where T is given by ( pTl| ) and (1.2), one can derive 
the following expression for the Kretschmann scalar 

1 

3" 



(22.1) K = Rap^sR"^^^ = CapjsC"^^^ + 2RapR''P - ^R" 



However, according to |18 , p. 19, we have 
(22.2) 



CaP'fSC"'^'^^ — 8{EapE"^ — HapH 



aP\ 
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where, relative to the frame Cq appearing in Lemma 21.1 , aU components of E and 
H involving cq are zero, and the ij components are given by 



H,,i = -'da' 



1 



where Sij is the same expression that appears in ( 21. (^ ), see p. 40 of Observe 
that in our situation, E and H are diagonal, since we are interested in the dcvel- 

" E^j/e^ and 



opments obtained in Lemma 21.2. It is natural to normalize Eij 
similarly for H. We will denote the diagonal components of Eij by E, 
to have expressions in S+, E_ and so on, and therefore we compute 

1 



We want 



1 



-JN2 



iV3)S- 

1 



E2 — E:i 



E2 + E-i 



N2{^+ + \/3S_) + ^{Ns - A^l)(S^ 



1 



3\/3 
2 



I]_(l - 2Sh 



2 



V3 



S-1 



S+(l 



9 



-N3){N2 



Observe that all other components of Ei and Hi can be computed from this, as Eij 
and Hij are both traceless. 

It is convenient to define the normalized Kretschmann scalar 
(22.3) k = Rc,p^sR°''^''^/0'^. 

The latter object can be expressed as a polynomial in the variables of Wainwright 
and Hsu. By the above observations and the fact that fl = 3fj./9^, we have 



i^[4+(37-2f]f72. 



We will associate a k and a RapR"^ to a solution to (2.1)-( pT3| ) in the following 
way. Since can be expressed in terms of the variables of Wainwright and 

Hsu, it is natural to define k by this expression multiplied by 0^, where obeys 
( pLl^ ). There is of course an ambiguity as to the initial value of 9, but we are only 
interested in the asymptotics, and any non-zero value will yield the same conclusion. 
We associate RapR"^ to a solution similarly. 

Lemma 22.1. The normalized Kretschmann scalar ( 22. j[ j is non-zero at the fixed 
points F, P^{II), at the non- special points on the Kasner circle, and at the type I 
stiff fluid points with 57 > 0. Consequently 

(22.4) limsup|K(T)| 



00 



for all solutions to (2.1)-(2.S) which have one such point as an a-limit point. 

Proof. The statement concerning the normalized Kretschmann scalar is a compu- 
tation. Equation (22.4) is a consequence of this computation, the fact that k = kO'^ 
and the fact that 6* — > 00 as r — > — cx), cf. (21.12). □ 

For some non- vacuum Taub type solutions with 2/3 < 7 < 2, the following lemma 
is needed. 
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Lemma 22.2. Consider a solution to with O > and 2/3 < 7 < 2 

such that 



(22.5) 

Then 



_lim (S+,I]_) = (-1,0). 



lim k(t) 



Proof. By Proposition p. 1| , the solution must satisfy S_ = and iV2 = N3. Observe 
that because of ( ^2.5 ), we have 17 — > 0, since ft decays exponentially for large, 

2. One can then prove that for 

e)r] 



cf. the proof of Lemma 



14.1 



Consequently, q 
any e > 0, there is a T such that 

(22.6) exp[(a^ + e)r] < rj(r) < exp[(a^ 

(22.7) exp[(6 + e)r] < Ni{t) < exp[(6 - e)r] 

(22.8) exp[(6 + e)r] < [iVi(iV2 + A^3)](r) < exp[(6 - e)r] 

(22.9) exp[(-6 + e)T] < 9^{t) < exp[(-6 - e)T] 

for all T <T, where = 3(2 — 7). However, the constraint can be written 
(1 - I]+)(l + E+) = n + - ^Ni{N2 + N3). 



By ( p2.6 )-(22.8), fl will dominate the right hand side, since it is non-zero. Since 
1 — E+ converges to 2, 1 + E+ will consequently have to be positive and of the order 
of magnitude fl. In particular, for every e > there is a T such that 



(22.10) 



exp[(a^ + e)T] < (1 + S+)(r) < exp[(a^ - e)r] 

q2 



S+)6' both diverge to infinity as r 



Observe that sin ce g-y < 4, flO^ and (1 
—00, by ( ^2.6| ), (22.9) and ( ^2.10 ). Other expressions of interest are Ni6'^ and 
Ni{N2 + N3)9'^ . The estimates ( [j2.6 )-(22.S) do not yield any conclusions concerning 
whether they are bounded or not. However, using (21.12), we have 



iVi(r)6|2(r) = Ari(O) 6*2(0) exp[/ (2 + q + AY.+)ds] 



Ni{0)e^{0)e^p[ (2(1 



-(37-2)f7 + 2S+(l + S+))ds] 



which is bounded since all the terms appearing in the integral are integrable by 
( ^2.6| ) and (22.10). A similar argument yields the same conclusion concerning 
Ni{N2 + N3)9\ 

Since the solution is of Taub type, we have Hi — NiY.^ and H2 — H3 — ~Hi/2. 
We also have E2 — E3 and 

2E2 = l^+a + S+) - ^N^ + ^NiiN2 + N3). 

Consequently the E field blows up and the H field remains bounded, and the lemma 
follows. □ 

Finally, we observe that RajsR"^ becomes unbounded in the matter case. 
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Lemma 22.3. Consider a solution to (2.1)-(2.S) with O > 0. Then 



lim R^bR"'^ 



Remark. How to associate R^pR"^ to a solution of (2.1)-(J2T3|) is clarified in the 



remarks preceding the statement of Lemma 22.1 
Proof. We have 



3p' = [1 + 3(7 - = -[1 + 3(7 - ir]n'e\ 



But by (|23|) and ( |21.12| ), we have 

n^iT)e^{T) = n^{o)e^{o) exp( / {-Aq + 2(37 - 2) + 4 + 4q)ds) 



{0)0^(0) exp(-37T), 



and the lemma follows. □ 



Lemma 22.4. Consider a class A development, not of type IX, with I = {t^,t-^) 
and 9 > Q. Then the corresponding solution to the equations of Wainwright and 
Hsu has existence interval M, and t ^ t± corresponds to t ^ ±00. 

Proof. The function 9 has to converge to infinity ast ^ t- for the following reason. 
Assume it does not. As 9 is monotone decreasing, we can assume it to be bounded 
on (t_,0]. By the constraint (21.8), aij and /i are then bounded on (i_,0], so that 
the same will be true of Uij by (21.4) and the fact that i_ > —00. But then one 
can extend the solution beyond contradicting the fact that / is the maximal 
existence interval. By ( |21.7| ), 6'-->0asi^oo = t+. Equation ( ^I.IC ) defines a 
diffeomorphism f : (t_,t+) — )■ (t_, t+), and w e get a solution to the equations of 
Wainwright and Hsu on (T_,r+). By ( 21.12 ), we conclude that the statement of 
the lemma holds. □ 

Lemma 22.5. Consider a Bianchi IX class A development with I (t_,t+) and 
1 < 7 < 2. According to Lemma 21.6, there is a to d I such that 9 > Q in I^ ~ 
(t^jto) and 9 < in 1+ = (to,i+)- The solution to the equations of Wainwright and 
Hsu corresponding to the interval /_ has existence interval (— cx3,r_), and t ^ t^ 
corresponds to t —00. Similarly, corresponds to (— oo,t+) with t ^ t^ 
corresponding to t ^ —00. 

Proof. Let us relate the different time coordinates on /_. According to equation 
( 21.1C| ), T has to satisfy dt/dr — 3/9. Define f (t) ~ J*^ 9{s)/Sds, where ti e 
Then t : I^ t(/_) is a diffeomorphism and strictly monotone on Since 9 is 
positive in f increases with t. 



Since 9 is continuous beyond to, it is clear that f (t) ^ r_ e M as t ^ io. To prove 
that t ^ t^ corresponds to r ^ —00, we make the following observation. One of 
the expressions 9 and d9/dt is unbounded on (i_,ti], since if both were bounded 
the same would be true of aij, p. and by ( ^1.7 ) and ( 21.4 ) respectively. Then 
we would be able to extend the solution beyond contradicting the fact that / 
is the maximal existence interval (observe that t- > —00 by Lemma 21.8). If f 
were bounded from below on /_, then 9 and 9' would be bounded on f((t_,ti]) 
by Lemma b3, and thus 9 and d9/dt would be bounded on (i_, ii]. Thus t ^ t- 
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corresponds to r — > — oo. Similar arguments yield the same conclusion concerning 
!+■ □ 
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